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Abstract 


This  thesis  describes  an  experimental  study  of  the 
electron-coupled  nuclear-nuclear  interactions  and  the 

Knight  Shift  in  white  tin  using  macroscopic  rotation 

* 

techniques.  A  description  is  given  of  a  high-speed  gas 

turbine  capable  of  operating  to  speeds  in  excess  of  5  kHz 

at  the 'magic'  angle  within  a  Varian  crossed-coil  probe. 

117  119 

Using  this  turbine  the  Sn  and  Sn  resonances 

were  observed  in  tin  metal.  The  lineshape  is  governed  by 

the  spin-lattice  relaxation  time  and  the  'pseudo-exchange' 

coupling  between  nuclei.  The  theoretical  lineshape  has 

been  synthesized  assuming  the  Ruderman-Kittel  model.  By 

comparing  the  experimental  and  theoretical  lineshapes  it 

was  found  that  the  nearest-neighbour  'pseudo-exchange' 

coupling  was  1.89  ±  0.09  kHz.  This  suggests  that  Alloul 

% 

and  Deltour  have  misinterpreted  their  data  in  obtaining  a 
value  of  4.1  kHz. 

The  room- temperature  isotropic  Knight  Shifts  of  the 

Sn  and  Sn^"^  isotopes  were  determined  to  be  0.724  +  0.003 

and  0.726  ±  0.003  respectively.  A  redetermination  of  the 

119  117 

ratio  of  the  Sn  and  Sn  magnetogyric  ratios  gave  the 
value  1.046535  +  0.000003. 
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CHAPTER  1.  INTRODUCTION 

When  an  ensemble  of  non-interacting  spin 
nuclei  is  placed  in  a  steady  magnetic  field  Hq,  the  two 
magnetic  energy  levels  of  the  system  are  separated  by 

AE  =  v  Ti  H  =  hv 
1  n  o 

where  y^  is  the  magnetogyric  ratio  of  the  nuclei .  The 
application  of  an  rf  magnetic  field,  of  frequency  v,  per¬ 
pendicular  to  the  main  field,  induces  transitions  between 
the  levels .  This  would  imply  that  the  nmr  spectrum  of  such 
a  system  consists  of  an  infinitely  sharp  line  at  frequency 
v.  If  the  nuclei  are  interacting  with  each  other,  the  nmr 
spectrum  ceases  to  be  a  sharp  line  and  becomes  spread  or 
broadened  over  a  band  of  frequencies . 

A  common  source  of  broadening  in  solids  arises 

from  the  nuclear  magnetic  dipole-dipole  interactions .  A 

nuclear  magnetic  dipole  produces ,  at  a  neighbouring  lattice 

site,  an  effective  additional  magnetic  field  with  an  angular 

2 

variation  of  (3  cos  0-1)  where  0  is  the  angle  between 
the  applied  field  and  the  internuclear  vector.  The 
presence  of  neighbouring  nuclei  can  thus  alter  the  resonance 
absorption  lineshape.  It  was  observed,  however,  that  in  many 
cases,  the  line  would  suddenly  narrow  when  the  sample  was 
heated.  This  was  taken  as  evidence  of  relative  motion 


* 
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between  nuclei  or  groups  of  nuclei  in  the  solid.  A 
given  neighbouring  nucleus  would  then  see  only  a  time 
averaged  field,  and  this,  in  general, is  smaller  than  the 
static  field.  Gutowsky  and  Pake  (1948)  showed  that  when 
the  frequency  of  motion  became  greater  than  the  width  of 
the  absorption  line  (measured  in  frequency  units) ,  narrowing 
occurred . 

Andrew,  Bradbury  and  Eades  (1958)  and  independently 
Lowe  (1959)  found  that  a  similar  narrowing  was  also  obtained 
when  the  sample  was  macros copically  rotated.  Further 
analysis  (Andrew  and  Newing,  1958)  showed  that  the  maximum 
narrowing  occurred  when  the  axis  of  rotation  was  at  the  so- 
called  ’magic'  angle  cos  ^(l//3)  to  the  main  magnetic  field. 

In  metals  there  is  an  additional  source  of  broaden¬ 
ing  above  that  of  the  dipolar.  This  arises  from  the  inter¬ 
action  of  the  nuclei  with  the  conduction  electrons.  This 
interaction  is  normally  stronger  than  the  dipolar  inter¬ 
action,  so  that  the  broadening  arising  from  it  dominates. 
There  are  three  types  of  broadening;  one  more  important  at 
low  temperatures ,  another  at  high  temperatures  and  the  third 
in  large  magnetic  fields.  The  presence  of  the  applied  mag¬ 
netic  field  polarises  the  conduction  electrons  in  the  metal. 
This  polarization  results  in  an  additional  magnetic  field 
at  the  nuclear  site,  giving  the  nucleus  a  shift  in  resonant 
frequency  relative  to  the  same  nucleus  in  a  non-metallic 
salt.  This  shift  is  known  as  the  Knight  Shift  (Knight, 
1949).  When  the  nucleus  is  ma  non-ciibic  environment,  this 


. 
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shift  depends  upon  the  angle  between  the  main  field  and 
the  crystallographic  axes  of  the  metal.  In  a  powder, 
therefore,  a  spread  of  resonant  frequencies  would  result 
from  the  crystallites  being  randomly  orientated.  This 
broadening  increases  in  proportion  to  the  applied  field. 

At  high  temperatures ,  the  electron-nuclear  interaction 
gives  a  finite  lifetime  to  the  nuclear  states .  This  results 
in  broadening  due  to  the  short  spin-lattice  relaxation 
time  T^.  As  the  relaxation  time  is  inversely  proportional 
to  temperature  (Korringa,  1950)  this  source  of  broadening 
is  relatively  unimportant  at  low  temperatures .  The  third 
effect  is  the  electron-coupled  nuclear-nuclear  interaction, 
which  at  low  temperatures  and  fields  has  a  dominant  effect 
on  the  metal  spectrum  (Bloembergen  and  Rowland,  1953; 
Rudermann  and  Kittel,  19  54)  .  This  arises  from  an  electron, 
scattered  by  its  hyperfine  interaction  with  one  nucleus, 
interacting  with  a  second  nucleus,  and  effectively  coupling 
the  nuclei  together.  Whilst  historically  this  effect  was 
first  seen  in  high  resolution  spectrometry  in  liquids 
(Hahn  and  Maxwell,  1952;  Gutowsky,  McCall  and  Slichter, 

1953)  it  has  since  been  noticed  in  many  metals.  In  general 
this  interaction  has  an  anisotropic  contribution. 

In  a  similar  fashion  to  the  narrowing  of  lines 
broadened  by  dipole  interactions,  the  macroscopic  rotation 
of  samples  leads  to  the  removal  of  other  anisotropic 
interactions.  This  enables  a  study  to  be  made  of  the  scalar 
part  (or  pseudo-exchange)  of  the  electron-coupled  nuclear- 


nuclear  interaction  without  the  necessity  of  considering 
a  spectrum  complicated  by  the  presence  of  any  anisotropic 
interaction . 

In  this  thesis  we  describe  an  experiment  conducted 
to  measure  the  scalar  'pseudo-exchange'  coupling  in  white 
(metallic)  tin.  Chapter  2  gives  details  of  the  theory  and 
experimentation  involved  in  the  method  of  macroscopic 
rotation.  Chapter  3  gives  a  brief  account  of  the  electron- 
coupled  nuclear-nuclear  interaction.  Chapter  4  gives  the 
analysis  of  the  data  obtained  from  tin.  The  results  are 
discussed  and  a  conclusion  given  in  Chapter  5. 
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CHAPTER  2.  MACROSCOPIC  ROTATION 
A)  Theory 

A  number  of  studies  have  been  undertaken  to 
examine  the  electron-coupled  interactions  in  tin  (steady- 
state  measurements:  Karimov  and  Shchegolev,  1961;  Sharma, 
Williams  and  Schone,  1969;  pulse  measurements:  McLachlan,  1968 
Alloul  and  Deltour,  19  69)  .  The  measurements  have  involved 
both  the  anisotropic  and  scalar  couplings.  Whilst  in 
principle  it  would  be  possible  to  remove  the  anisotropic 
interactions  with  a  series  of  pulses  (cf  removal  of 
dipolar  line  broadening:  Mansfield  and  Ware,  1966;  Ostroff 
and  Waugh,  1966)  ,  in  practice  there  are  difficulties. 

If  measurements  were  made  using  a  single  metal  crystal 
there  would  be  the  difficulty  of  overcoming  signal  distor¬ 
tions  due  to  the  finite  skin  depth.  On  the  other  hand,  if 
a  powder  was  used,  then  the  line  shape  would  be  complicated 
by  the  presence  of  the  anisotropic  Knight  Shift.  This  is 
a  macroscopic  rather  than  a  microscopic  broadening  of  the 
line  and  therefore  cannot  be  removed  by  a  system  of  pulses 
(Andrew,  1970) .  Removing  the  anisotropic  broadening  by 
rotation  therefore  offers  certain  advantages. 

It  is  interesting  to  note  that  historically  the 
first  measurements  involving  sample  rotation  were  not 
made  for  the  specific  purpose  of  narrowing  the  line. 

Instead,  Andrew,  Bradbury  and  Eades  (1958)  used  the 


6 


technique  to  confirm  the  theoretical  predictions  of 
Anderson  (1954)  and  Pake  (1956)  that  internal  motion  in  a 
solid  narrowed  the  dipolar-broadened  spectrum  in  such  a 
manner  that  its  second  moment  remained  constant,  whilst 
its  fourth  moment  increased.  When  Andrew  et  al.  rotated 
a  sample  of  NaCl,  they  found  that  the  central  part  of  the 
spectrum  narrowed  but  satellite  lines  appeared.  These 
were  separated  from  the  main  line  by  multiples  of  the 
rotation  frequency.  This  increase  in  intensity  in  the  wings 
keeps  the  second  moment  invariant  whilst  increasing  the 
fourth  moment. 

The  effect  of  rotation  is  best  seen  by  considering 
the  simple  case  of  an  ensemble  of  pairs  of  nuclei  with 
spin  —  arranged  in  a  static  array  with  their  internuclear 
vectors  parallel  to  each  other.  The  nmr  spectrum  of  each 
pair  of  non-rotating  nuclei  consists  of  two  lines  of  equal 
intensity  at  the  angular  frequencies 

=  Ya [Hq  +  ey^R  3^3  cos20  “1)3/  (D 

where  y  is  the  magnetogyric  ratio  of  the  nuclei  at 

resonance,  Hq  is  the  constant  applied  field,  y^  is  the 

magnetic  moment  of  the  nuclear  neighbour  in  each  pair,  R 

is  the  internuclear  vector,  0  is  the  angle  between  H 

-  o 

and  R  and  e  is  3/2  if  the  nuclei  in  the  pair  are  identical 
and  1  if  they  are  not  (Pake,  1948) .  Thus  the  spectrum 
of  this  static  array  can  be  considered  to  arise  from  the 
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presence,  at  the  site  of  the  nuclei  at  resonance,  of  an 
effective  local  field  whose  component  along  Ho  is  given  by 

Hloc  =  —  cos^S  "  1)  (2) 

-3 

where  c  =  e  R  (3) 

Instead  of  considering  a  static  array,  let  the  nuclei 
be  capable  of  being  rotated  about  an  axis  making  an  angle 
a  to  H^,  and,  for  simplicity,  perpendicular  to  the  inter- 
nuclear  vector  of  the  pair  (see  figure  (2.1a)).  Imagine 
that  when  there  is  no  rotation,  the  nuclear  pairs  are 
distributed  with  random,  but  fixed,  azimuth  angle  <j>  about 
the  rotation  axis.  The  spectrum  for  each  pair,  with  a 
given  <j> ,  consists  of  two  lines  given  by  equation  (1)  ,  in 
which  we  set 


cos  0  =  sin  a  cos<J)  .  (4) 

The  spectrum  for  all  pairs  ,  taking  all  4>  with  equal 
probability,  has  a  second  moment 


S„  =  —  f  c^v  ^(9  sin^acos%  -  6  sin^acos^cj)  +  1)  d<J> 
z  2  ti  a 

o 


2  2,27  .4  ,  .  2  , 

c  y  (-3—  sin  a  -  3sm  a  +  1)  . 
a  o 


(5) 


If  all  the  pairs  are  now  allowed  to  rotate  with  the  same 
uniform  angular  velocity  co^,  we  then  have 


<f>  +  03  t 

o  r 


(6) 


. 
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Figure  (2.1) 


(a)  A  system  of  two  nuclei  A  and  B  capable  of  being  rotated 

about  an  axis  at  an  angle  a  to  the  main  magnetic  field  H^. 


(b)  Nmr  spectrum  for  the  above  pair  of  nuclei  when  they  are 

rotating  at  a  rate  to  about  an  axis  perpendicular  to  their 


internuclear  vector.  The  frequencies  ojq  and 

1  2 

and  2  c  Ya  (3  cos  a  -  1)  respectively. 


Y  H 
'  a  o 


are 


3 


J^(m) 

_j _ i _ 

(X>a±00,-4-^r 


Jo(m)  ( b) 


jgCm) 

1 _ L. . . 

C0o±  cO,  -2oOr  cOo  ±  60 1  ~  ^1  +2- 00 r  <4>  -  0) 

ANGULAR  FREQUENCY 
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so  that  each  nucleus  experiences  a  local  field  which 
consists  of  static  and  fluctuating  components.  The  Larmor 
frequency  of  the  nuclei  is  then  given  by  (Andrew  and 
Newing,  1958) 

w  =  Y  (H  +  H,  ) 

'a  o  loc 


where 


=  to  +  w,  +  w~  cos  2  to  t 
o  1  —  2  r 


(7) 


w 

o 


H 

o 


to 


w 


1 

2 


1  n  2 

•?r  c  y  (3  cos  a  -  1) 

2  1  a 


3  .  2 

7T  c  y  sin  a 

2  1  a 


(8) 


We  note  that  the  Larmor  frequency  is  frequency -modulated 
at  2tor«  Thus  the  spectrum  consists  of  two  carrier  com¬ 
ponents  at  frequencies  toQ  +  to^ ,  each  of  which  has  an  infi¬ 
nite  set  of  sidebands  situated  at  multiples  of  2tor  (see 

figure  (2.1b)).  The  intensity  of  the  nth  sideband  relative 

2  2 

to  that  of  the  central  component  is  Jn  (m)/JQ  (m)  where 
Jn (m)  is  the  Bessel  function  of  order  n,  and  m  is  given  by 

m  =  w2/2tor.  (9) 

The  second  moment  of  this  spectrum  is 

+00  9  2 

=  =■  E{(w1  +  2noj  )  J  (m) 

2  2  1  r  n 

n=— co 

+  (-to,  +  2noo  2)  ^  J  2  (m)  } 

1  r  n 


(10) 


*- 
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which,  by  making  use  of  the  theorems  of  Watson  (1922), 
becomes 


Making  use  of  equations  (8)  and  (11)  one  finds  that  the 

second  moment  is  identical  to  that  of  the  static  array. 

As  to  increases ,  the  sidebands  move  out  and  become  less 

intense.  If  we  ignore  their  contributions,  the  spectrum 

then  consists  of  two  carrier  frequencies.  The  second 

2 

moment  of  these  is  to-^  ,  which  is  smaller  than  the  static 
second  moment.  We  can  generalize  this  analysis  to  a  more 
random  array  of  nuclei,  whose  spectrum  is  broad  rather 
than  consisting  of  delta  functions,  and  see  that  rotation 
would  lead  to  a  narrower  line.  We  note  from  equation  (8) 
that  if 

a  =  cos""1  (1//3 )  =  54  °  44  1  (12) 

the  centre  lines  collapse  into  a  single  line  at  oj  ,  and 
the  dipolar  second  moment  (excluding  the  contribution  of 
the  sidebands)  reduces  to  zero. 

We  shall  now  proceed  to  consider  macroscopic  rota¬ 
tion  in  a  more  general  fashion.  Theoretical  considerations 
have  been  made  on  the  narrowing  of  dipolar-broadened  lines 
(Andrew  and  Newing,  1958;  Andrew  and  Jenks,  1962)  and  on 
the  removal  of  the  anisotropic  Knight  Shift  (Schwind,  1967 
but  the  most  general  comprehensive  study  has  been  undertaken 


' 
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by  Andrew  and  Farnell  (1968)  who  were  able  to  point  out  a 
limitation  to  the  rotation  technique  in  removing  aniso¬ 
tropic  effects . 

The  truncated  Hamiltonian  which  determines  the  nmr 
spectrum  is  ,  excluding  quadrupole  interactions, 


H  =  H  +  [  H]  +Ht. 

z  D  trunc.  J 


(13) 


H  is  the  Zeeman  interaction  between  the  nuclei  and  the 
z 

applied  field: 


H  =  -  £  Y .  I  . 
z  1  1  —1 

l 


(1  -  K. )  .  H 

=  =1  — o 


(14) 


where  I_^  is  the  spin  operator  of  the  nucleus  i,  y^  is  the 
magnetogyric  ratio,  and  is  the  Knight  Shift  tensor.  HD 
is  the  truncated  dipolar  interaction  (Van  Vleck,  1948; 
Slichter ,  1963)  : 


[ H  ]  =  E  h  I  .  .  D,  .  I , 

D  trunc.  .  .  —l  =t  — d 

i<3 


(15) 


where  is  the  truncated  dipolar  interaction  tensor,  all 

2 

terms  of  which  contain  the  factor  (3  cos  0. .  -  1) ,  with 
0^j  being  the  angle  between  Hq  and  the  internuclear  vector 


R.  . 
“13 


H  is  the  indirect  electron-coupled  nuclear-nuclear 
J 


spin  interaction: 


Ht  =  E  hi.  .  J.  .  .  I  .  (16) 

J  .  .  —i  =n 

i<3 

where  J. .  is  the  indirect  interaction  tensor. 

=  13 

When  the  sample  is  rotated,  the  Hamiltonian  becomes 
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time-dependent.  It  can  then  be  expressed  as 

H  =  H  +  H  ( t)  (17) 

where  H  is  the  time-averaged  Hamiltonian,  and  H(t)  is  the 
time-dependent  interaction  whose  average  value  is  zero. 

H(t) ,  which  is  small  compared  to  H  when  the  main  field  is 
strong,  is  responsible  for  the  generation  of  the  satellite 
lines  discussed  above.  H,  which  determines  the  central 
part  of  the  spectrum,  is  given  by  the  sum  of  H  ,  E^  and 
H  .  We  shall  consider  the  case  where  the  satellites  are 
sufficiently  far  removed  from  the  centre  of  the  spectrum, 
so  that  they  are  not  observed;  i.e.  only  the  time-averaged 
Hamiltonian  will  be  examined.  Since  the  Knight  Shift 
components  are  so  small,  only  the  components  along  Hq  need 
to  be  considered  (Abragam,  1961)  and  H  may  be  written  as 


H  =  -  ft  E 


Y±(l  - 


K.  ) 

1ZZ 


1Z 


H 

o 


+ 


E 

a3 


E  (h 

i<j 


J..  0  +  D .  .  0 ) I .  I.0 

ija3  ija3  la  j3 


where  a  and  3  are  suffixes  for  cartesian  axes 
frame  and 


(18) 

in  the  lab 


K  + 


K  ' 

z 


+  K 


x 


K  ' 

y  2 


(19) 


with  X  being  the  direction  cosines  of  the  principal  axes 

P 

of  the  anisotropic  Knight  Shift  relative  to  H^,  K^  is  the 
pth  component  of  the  anisotropic  Knight  Shift  tensor  K' 
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and  K  is  the  isotropic  Knight  Shift. 

Now,  when  the  sample  is  rotated  with  an  angular 
velocity  about  an  axis  inclined  at  an  angle  to 
and  at  an  angle  y^,  and  y^  to  the  principal  axes  of  K '  , 
we  have  that 


X  .  i =  cosa.cosy  . , 

1  O  1 


sincusiny  .icos  (co  t 

O  3.  3T 


(20) 


where  tJk,  is  the  azimuth  angle  of  the  ith  principal  axis  at 
t  =  0.  Using  eqn.  (20)  in  equation  (19)  we  find 


K  =  K  +  ^  sin2a0 (K '  +  K'  +  k'  ) 

zz  2  3  x  y  z 


+  ~(3cos2a0  -  1)  E  K.',  cos2i p./. 

m  •  f  1  jL 

1 


(21) 


When  a ^  has  the  special  value,  cos  ^(l//3)  ,  the  so-called 
magic  angle,  this  mean  value  reduces  to 


K  =  K  (22) 

zz 


asK'  +  K '  +  K '  =0. 

x  y  z 

The  dipolar  interaction  is  treated  in  a  similar 

fashion ; 


E 

i>  j 


E  D.  .  0  I  .  I  .  0 
ae  13“e  101  36 


,  0  2  3n-1/0  2  Q  1N  (I.. I. 

=  E  y  .  y  .  (  8tt  r.  .  )  (3cos  0..-1)  — l  —  3 

i>j  1  3  1D  J 


-  31.  I 


1Z 


(23) 


Since 


cos  0 .  .  =  cosa0cosfi.  . 
13  3  13 


sina^sinQ.  .cos  (to  t  +  d) .  . ) 
3  13  r  13 


(24) 
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Figure  (2.2) 


A  diagram  of  a  general  system  of  rotating  axes  to  illustrate  the 
nomenclature  used  in  Chapter  2. 
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where  ^  is  the  angle  between  the  axes  of  rotation  and 

R.  .  ,  with  (J> .  .  being  the  azimuth  angle,  at  zero  time,  of 
13  1 D 

R.  .  ,  we  find  that 
-ID 


2  12  2  1 
cos  0.  .  =  -^(3cos  cu  -  1)  (3cos  ft.  .  -  1)  +  — . 
ID  o  3  13  3 


(25) 


Thus?  at  the  'magic  angle', 

(Scos2©^  -  1)  =  0  (26) 

and  the  dipolar  contribution  vanishes . 

The  J  tensor  must  be  treated  differently.  It  may 
be  split  into  three  parts ; 

J  =  Jl  +  J*  +  J1"  (27) 

where 

J  =  1/3  trace  J 
g*  =  \  U  +  g')  -  Ji 
j+  =  |  (j  -  j') 

in  which  J'  is  the  transpose  of  J.  Clearly  this  divides 

J  into  an  isotropic  part  Jl,  a  symmetric  part  J*  and  an 

t 

antisymmetric  part  g  . 

In  tin,  symmetry  considerations  show  that 
t 

J  should  be  zero  (see  Appendix  II)  .  (in  fact;  no  case 

has  yet  been  observed  in  solids  where  this  is  not  the  case 

*1" 

(Andrew,  1970)  .)  We  shall  therefore  take  J  as  zero. 


2 


. 


The  component  J1  is  unaffected  by  the  rotation.  The  term 


can  be  written  as 


E 

a$ 


i  ja£ 


I. 

la 


E  E 
n  a3 


C.  .  C.  . 
ljan  13 


Q  l.  i  .  DJ.  . 
3n  1a  33  ljn 


(28) 


where  C  are  the  direction  cosines  of  the  nth  principal 
axis  of  the  tensor  relative  to  the  laboratory  axes. 

For  each  set  (i,  j)  there  are  27  terms,  not  all  of  which 
correspond  to  transitions  near  the  main  resonance  line. 
Truncating  as  for  the  dipolar  term  we  obtain 


J?.  QI. 
liaB  1 


+ 


)  (C 


2 

i  j  32 


i  j  31 


I 

z  3 


)  . 


When  the  crystal  is  rotated,  the  direction  cosines  become 
time-dependent  and  the  average  value  of  their  squares  is 


C?.0  =  i  (3cos2a0-l)  (3cos2y.  .  -1)  +  i 

i33n  6  3  Ki^n  3 


(30) 


where  u . .  are  the  angles  between  the  axis  of  rotation  and 
ljn  r 

the  principal  axes  of  the  tensor  . .  Thus,  at  the  magic 
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angle 


[  Z  E  J 
i< j  a3 


(31) 


Collecting  terms,  we  see  that  on  macroscopic  rota¬ 


tion  at  the  magic  angle  the  resonance  line  is  determined 
by  the  time-averaged  Hamiltonian 


H  =  -  E  y.  fi  I  (l-K.)H  +  I  J .  .  h  I .  .  I  .  . 

i  i  zi  i  o  i<j  13  -1  -3 


(32) 


We  find  that  the  nmr  spectrum  is  determined  by  the  Knight 
Shift  K  and  the  coupling  constant  J  as  in  the  case  of  an 
isotropic  liquid.  (In  a  liquid  the  Knight  Shift  is 
replaced  by  the  scalar  chemical  shift.)  We  shall  attempt 
to  analyse  the  data  using  a  Hamiltonian  of  this  form  in 
chapter  4. 

Andrew,  Carolan  and  Randall  (1971a,  1^71b)  have 

undertaken  the  only  other  study  of  this  kind.  Their 

63  6  5 

measurements  were  made  on  Cu  and  Cu  in  metallic  copper 
and  the  value  of  the  J  coupling  was  obtained  by  measuring  the 

second  moment  of  the  resonance  absorption.  However,  as  will  be 
seen  in  Chapter  4}  we  have  used  an  entirely  different 
approach  in  obtaining  the  value  of  the  J  coupling. 
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B)  Experimentation 

To  obtain  a  spectrum  determined  by  the  'liquid- 
like'  Hamiltonian,  it  is  necessary  to  rotate  the  sample 
relative  to  the  external  field.  This  may  be  achieved  by 
rotating  either  the  sample  or  the  field.  The  field  may  be 
rotated  electrically  much  faster  than  the  sample  can  be 
mechanically  spun.  Bradbury,  Eades  and  McCarten  (1968) 
report  generating  a  270  gauss  field  rotating  at  27kHz. 
However,  the  problem  here  is  to  create  a  field  that  is 
large  enough  to  be  useful.  Overcoming  this  problem,  by 
spinning  the  sample,  limits  the  speeds  to  about  10kHz.  The 
sample  is  normally  held  in  some  type  of  container  which  is 
suspended  upon  a  set  of  air  bearings.  This  means  that 
if  the  sample  is  slightly  unbalanced,  it  can  choose  its 
own  axis  of  rotation  and  spin  freely . 

If  the  sample  was  mounted  on  a  rigid  shaft  as  is 
shown  in  figure  (2.3a)  then  a  slight  imbalance  would,  at 
high  speeds ,  create  such  large  forces  that  the  mechanical 
bearings  supporting  the  shaft  would  quickly  fail.  Instead 
of  making  use  of  air  bearings^,  it  is  possible  to  rotate  a 
sample  that  is  suspended  on  magnetic  bearings  (Beams,  1964). 
However,  this  type  of  suspension  is  not  suitable  because, 
even  if  the  sample  was  sufficiently  magnetic  to  be 
supported,  the  need  to  use  a  magnetic  field  gradient  is  in 
conflict  with  the  homogeneous  field  required  in  nmr;  the 
axis  of  rotation  is  normally  in  such  a  direction  that 
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Figure  (2.3) 


(a)  Schematic  of  stiff-shaft  turbine. 


(b)  Simplified  version  of  air  turbine. 


The  thrust  gases  are  issuing  at  a  velocity  v. 
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the  field  gradients  would  not  be  averaged  out.  Furthermore, 
this  fact  also  implies  that  the  rotation  would  not  be  at 
the  'magic  angle'. 

Details  of  an  early  air  turbine  were  given  by 
Henriot  and  Huguenard  (1925)  .  A  simple  adaptation  of  their 
turbine  (Beams,  1937)  is  shown  in  figure  (2.3b) .  The  tur¬ 
bine  or  rotor  is  supported  by  a  series  of  air  jets  in  the 
lower  support  or  stator.  These  jets  are  angled  away  from 
the  vertical  in  such  a  manner  as  to  provide  both  support 
and  propulsion.  The  turbine,  when  spinning  steadily,  is 
very  stable  and  will  actually  operate  upside-down.  However, 
when  it  starts  it  does  not  normally  spin  centrally  on  the 
stator  but  instead  tends  to  bounce  around.  This  motion 
must  be  damped  out  before  the  turbine  will  achieve  its 
highest  speeds.  The  most  efficient  and  easiest  method  of 
doing  this  is  to  use  the  tips  of  the  fingers.  Mechanical 
self-starting  was  tried,  but  when  it  was  used  needed 
constant  readjustment  due  to  turbine  wear.  (Fingers  have 
the  advantage  of  being  naturally  self-repairing.) 

The  turbine  system  used  in  the  measurements  is  an 
adaptation  of  that  used  by  Beams  (1937)  to  support  a  rotating 
mirror.  It  was  modified  to  allow  its  use  inside  the  Varian 
crossed-coil  probes  provided  with  the  Varian  nmr  spectro¬ 
meter.  It  was  felt  that  as  the  sample  need  not  rotate  in 
the  vertical  position,  it  would  be  preferable  to  orient  it 
at  54° 44'  to  a  horizontal  magnetic  field  rather  than  risk 
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straining  the  magnet  yoke  in  orientating  the  field  at  the 
magic  angle  to  a  vertically  spinning  sample  (contrast 
Andrew,  1970) .  The  fact  that  the  sample  spins  stably  can 
be  explained  in  terms  of  Bernoulli's  principle.  Referring 
to  the  simplified  version  of  the  air  turbine  and  bearing 
shown  in  figure  (2.3b)  we  see  that 


jp  +  pv  /2  =  P 
lower  surface  ' 


atmos 


(20) 


and  therefore  the  atmospheric  pressure  prevents 
the  turbine  from  being  ejected  from  the  stator  by 
the  propelling  gas.  If  we  now  suppose  that  the  turbine 
shifts  away  from  the  centre  of  the  stator,  then  there  is 
a  decrease  in  the  velocity  of  air  flowing  through  the  gap 
between  the  turbine  and  stator.  The  resulting  pressure 
increase  recentralizes  the  rotor. 

Details  of  the  construction  of  the  rotor  are  shown 
in  figure  (2.4)  .  The  turbine  (figure  (2.4a))  is  constructed 
from  nylon.  Great  care  and  precision  is  needed  to  insure 
that  the  flutes  are  all  milled  equally  deep,  uniformly  wide 
and  correctly  placed  on  the  surface  of  the  rotor  relative 
to  the  air  holes  in  the  stator.  (Despite  this  careful 
construction,  there  was  often  a  failure  rate  of  one  in  two.) 
The  hollow  in  the  centre  of  the  turbine  is  used  to  contain 
the  sample.  It  is  necessary  to  construct  the  stator  in 
two  parts.  The  upper  piece  (figure  (2.4c))  contains  the 
jets  to  spin  the  turbine. 
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Figure  (2.4) 

Construction  of  the  turbine  used  in  the  macroscopic  rotation 
of  the  nmr  sample. 

(a)  Diagram  of  the  turbine. 

(b)  Diagram  of  the  lower  part  of  the  stator. 

(c)  Diagram  of  the  upper  part  of  the  stator  showing  position 
of  air  jets . 

(d)  Diagram  showing  the  total  turbine  assembly  as  used 
experimentally . 
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The  piece  shown  in  figure  (2.4c)  is  set  on  top  of 
the  lower  stator  part  (figure  2.4b))  at  an  angle  of  about 
50°  to  the  central  axis  of  the  lower  piece.  Rotating  the 
whole  stator;  therefore,  allows  the  axis  of  rotation  of  the 
turbine  to  be  turned  so  that  it  makes  an  angle  between  50° 
to  90°  to  a  magnetic  field  applied  perpendicular  to  the  sta¬ 
tor  axis.  The  two  stator  pieces,  which  are  made  of  nylon, 
are  fused  together,  after  cleaning  with  trichloroethylene, 
with  formic  acid.  Provided  that  the  two  surfaces  have  not 
been  machined  'optically  flat' ,  this  creates  a  very  effective 
seal.  After  fusing,  the  pieces  are  machined  to  fit  down  the 
3/4"  sample-hole  of  the  Varian  probe.  The  assembled  stator 
is  mounted  on  a  piece  of  perspex  tubing  as  is  shown  in 
figure  (2. 4d) .  This  is  necessary  as  it  is  not  normally 
possible  to  start  the  turbine  inside  the  probe.  It  is 
therefore  started  outside  the  probe  and  then  brought  down 
level  with  the  receiver  coils  inside  the  probe.  Unfortunate¬ 
ly,  this  type  of  motion  is  one  against  which  the  turbine 
is  the  least  stable.  If  the  movement  is  jerky,  the  rotor 
will  touch  the  stator  and  fly  out.  As  the  turbine  circum¬ 
ference  is  often  travelling  near  the  velocity  of  sound  in 
the  propelling  medium,  the  turbine  travels  a  great  distance 
along  the  floor.  It  is  therefore  advisable  to  enclose  the 
magnet  to  prevent  loss  of  the  rotor.  The  enclosure  can 
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also  provide  increased  thermal  stability  for  the  magnet. 

Of  the  turbines  that  spin,  only  50%  will  do  so  inside  the 
probe.  This  is  apparently  due  to  the  unbalancing  effect  of 
the  air  stream  from  the  stator  being  reflected  back  onto 
the  turbine  from  the  probe  walls . 

The  maximum  theoretical  rotation  rate  for  the  tur¬ 
bine  occurs  when  its  circumference  is  travelling  at  the 
velocity  of  sound  in  the  propelling  medium.  Normally  this 
maximum  is  not  reached  as  a  number  of  speed-limiting  factors 
are  present.  Imbalance  will  cause  the  rotor  to  spin  un¬ 
stably  or  'jump'.  Frictional  forces  acting  on  the  sides  of 
the  rotor  will  counterbalance  the  thrust  provided  by  the 
jets.  The  stator  support  may  vibrate  causing  the  turbine 
and  stator  to  touch.  In  addition,  when  the  propulsion 
pressure  is  high,  it  is  likely  that  there  will  be  a  certain 
amount  of  turbulence  in  the  pipes  supplying  air  to  the 
turbine,  which  is  capable  of  upsetting  the  turbine.  It 
should  therefore  be  expected  that  Helium  gas  would  generate 
higher  speeds  than  Nitrogen  or  Air.  This  was  found  to  be 
the  case  experimentally.  Figure  (2.5)  shows  the  variation 
of  rotation  speeds  measured  as  a  function  of  thrusting  press¬ 
ure  for  both  Helium  and  Nitrogen.  Hydrogen  is  even  more 
effective  but  because  of  the  danger  involved,  it  is  not 
normally  used.  (Sparks  may  be  generated  when  the  stator 
and  turbine  touch.) 

The  actual  rotation  rates  are  easy  to  measure. 
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Figure  (2.5) 

Variation  of  rotational  rate  as  a  function  of  thrust  when 
Helium  and  Nitrogen  are  used  as  propelling  gases. 
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At  low  speeds,  it  is  possible  to  use  a  stroboscope  to  'stop' 
the  motion.  At  speeds  much  greater  than  100,000  rpm  it  is 
not  easy  to  determine  at  what  harmonic  of  the  stroboscope 
flash  rate  the  turbine  is  turning.  If  the  turbine  is 
painted  half  matt-black  and  half  silver,  then  light  reflected 
from  the  turbine  onto  a  photo-trans is  tor  creates  sufficient 
signal  to  trigger  a  counter.  However,  when  the  turbine  is 
sitting  at  an  angle  deep  inside  the  probe,  there  is  in¬ 
sufficient  light  reflected  back  to  the  transistor.  On 
hearing  the  turbine  rotate,  it  becomes  obvious  what  is 
another  method  of  measuring  the  rotation  rate.  Small 
imperfections  on  the  rotor  surface  set  up  an  audio  signal 
with  tones  equal  to  the  fundamental  rotation  rate  and  its 
harmonics.  (We  have  detected  up  to  the  12th  harmonic.) 

At  high  rotation  rates  the  harmonics  are  not  perceived  and 
it  is  then  a  simple  manner  to  beat  the  signal  against  that 
of  an  oscillator  which  feeds  both  a  counter  and  a  loud¬ 
speaker.  At  very  high  speeds  even  the  fundamental  may  be 
difficult  to  detect.  The  greatest  speed  obtained  was 
12kHz  (720,000  rpm) ,  but  this  was  with  a  solid  rotor  rather 
than  one  containing  a  sample.  That  particular  turbine 
'blew'  up  at  a  later  time  when  'turning  over'  at  7.2kHz. 

If  the  sample  is  nonmetallic,  it  may  be  packed 
tightly  into  the  centre  of  the  turbine  and  capped.  At  low 
speeds,  plastic  wood  is  very  useful  but  at  high  speeds 
the  edges  of  the  turbine  begin  to  creep,  allowing  the  cap 
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to  shift.  This  either  unbalances  the  turbine  and/or 
scatters  the  sample.  A  better  method  is  to  mix  the  sample 
with  an  epoxy  which  is  then  allowed  to  set  in  the  turbine. 

A  good  room- temperature-setting  epoxy  is  STY CAST  2850FT 
manufactured  by  Emerson  and  Cuming,  Inc.,  Canton,  Mass., 
U.S.A.  When  used  in  conjunction  with  catalyst  24,  it  is 
very  liquid  with  a  pot  life  of  about  45  mins.  The  fact 
that  the  epoxy  flows  freely  means  that  a  large  amount  of 
sample  may  be  mixed  in  before  the  epoxy  is  too  stiff  to 
pour.  To  prevent  the  formation  of  air  bubbles,  which 
would  unbalance  the  turbine,  it  is  placed,  after 
filling^in  a  bell  jar,  which  is  then  evacuated.  This  improves 
the  number  of  turbines  that  rotate  after  filling  to  50% from 
the  10%  that  would  spin  before  the  evacuation  technique 
was  used.  This  number  could  be  even  further  increased  if  - 
less  sample  was  mixed  in  with  the  epoxy.  The  decrease  in 
sample  volume  is,  however,  undesirable,  as  it  decreases  the 
filling  factor  and  hence  the  s ignal-to-noise  ratio.  If  an 
enriched  sample  was  to  be  mixed  in  with  the  epoxy,  it  would 
be  advisable  to  leave  the  epoxy  very  liquid  to  avoid  wasting 
the  sample.  Using  epoxy  rather  than  plastic  wood  offers 
a  large  number  of  advantages .  The  epoxy  has  very  good 
mechanical  strength.  Even  if  the  turbine  breaks  up,  the 
epoxy  will  hold  together,  preventing  the  sample  from  con¬ 
taminating  the  probe.  When  a  metal  powder  sample  is  used, 
the  high  resistance  of  the  epoxy  efficiently  isolates  the 
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particles  from  each  other,  alleviating  any  problem  of  skin 
depth.  The  epoxy  is  distributed  evenly  throughout  the 
turbine  rather  than  sitting  on  top  as  a  cap  would.  This, 
together  with  the  fact  that  the  epoxy  binds  onto  the  edges 
of  the  nylon  preventing  them  from  spreading,  leads  to  a 
greater  mechanical  stability  at  high  speeds. 

The  problem  of  a  poor  filling- factor ,  when  using 
a  sample-epoxy  mix,  could  be  avoided  if  it  were  possible  to 
machine  the  whole  turbine  from  the  mix.  This  was  tried,  but 
it  was  found  that  the  epoxy,  especially  when  metal  particles 
were  added,  was  so  abrasive  that  it  was  impossible  to  machine 
the  turbines  to  any  accuracy.  To  overcome  the  signal- to- 
noise  problem,  it  was  necessary  to  employ  signal  averaging. 

The  actual  experimental  set-up  is  shown  as  a  block 
diagram  in  figure  (2.6).  The  turbine  is  driven  with  air 
supplied  by  a  compressor  with  an  operating  output  pressure 
fluctuating  between  90  and  110  Ibs/sq  in.  By  the  use  of  a 
series  of  regulators  and  storage  tanks,  the  air  flow  to  the 
turbine  could  be  regulated.  Over  a  day  the  rotation  rate 

4 

was  stable  to  around  1  part  in  10  .  To  gain  the  stability 
needed  for  long-term  signal-averaging,  the  frequency  of  the 
Varian  spectrometer,  operating  at  6MHz,  is  phase-locked 
to  a  1MHz  frequency  standard  and  is  stable  to 
1  part  in  10  .  The  magnet  field  is  controlled  by  the 
Varian  Mark  II  Fieldial  via  a  Hall  probe  sensing  element. 

We  have  attached  a  blackened  disk  with  small  silver  markings 
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Figure  (2.6) 

Block  diagram  of  experimental  arrangement  used  when  nmr 
samples  are  being  macros copically  rotated. 
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to  the  shaft  of  the  Fieldial  potentiometer  used  to  sweep 
the  field.  By  the  reflection  of  light  from  these  markings 
onto  a  phototransistor,  a  trigger  signal  is  generated  for 
controlling  the  Fabri-Tek  1062  instrument  computer.  The 
generated  signal  was  squared  up  to  give  a  pulse  with  a  good 
rise-time  so  that  the  computer  would  trigger  consistently 
at  the  same  point  in  the  field  sweep.  The  microswitch 
provided  by  Varian  in  the  Fieldial  was  found  to  be  insuf¬ 
ficiently  accurate  for  this  purpose.  Pulses  from  points  on 
the  disk  other  than  the  starting  position  could  be  fed  into 
the  computer  and  superimposed  on  the  resonances  to  provide 
calibration  markings.  The  sweep  rates  of  the  magnet  field 
were  altered  to  match  up  with  the  available  sweep  times  on 
the  Fabri-Tek.  This  provides  for  the  most  efficient  use  of 
the  computer;  its  ‘dead  time’  relative  to  the  field  sweep  time 
is  less  than  6%.  The  nmr  signal  is  audio-detected  in  the 
Varian  spectrometer  and  then  fed  into  a  narrow-band 
amplifier  and  lockin  detector  of  the  type  described  by 
Schuster  (1951)  . 

In  the  case  of  tin,  it  was  often  necessary  to  make 
2048  runs  over  a  period  of  48  hrs  to  achieve  an  adequate 
signal- to-noise  ratio. 


31 


CHAPTER  3.  THEORY  OF  THE  ELECTRON-COUPLED 
NUCLEAR-NUCLEAR  INTERACTION 

In  a  study  made  on  thallium  oxide ,  Bloembergen  and 

Rowland  (1953)  found  that  the  linewidth  of  one  of  the 

thallium  isotopes  was  much  greater  than  that  of  the  other. 

203  205 

The  extra  width  of  T1  compared  to  T1  was 

surprising  as  their  magnetic  moments  only  differ  by  1%/  so 

that  it  could  be  expected  that  both  isotopes  would  interact 

with  their  surroundings  to  the  same  extent.  The  only  manner 

in  which  the  isotopes  differ  greatly  is  in  their  natural 

20  3 

abundances.  The  T1  is  30%  abundant  compared  to  70%  for 
205  203 

the  Tl  ,  so  that  a  Tl  nucleus  will  on  average  have  a 

20  5 

larger  number  of  unlike  neighbours  compared  to  a  Tl 
nucleus .  The  experimental  data  may  be  explained  if  there 
is  an  exchange  interaction.  Exchange  interactions  between 
unlike  nuclei  broaden  the  resonance  line  whilst  those 
between  like  neighbours  narrow  it  (Van  Vleck,  1948) .  Hahn 
and  Maxwell  (1952)  and  independently  Gutowsky,  McCall  and 
Slichter  (1953)  have  shown  that  there  exists  a  coupling 
between  the  magnetic  moments  of  two  nuclei  in  a  molecule, 
via  their  hyperfine  interaction  with  the  surrounding 
electrons.  Ruderman  and  Kittel  (1954)  showed  that  a  simi¬ 
lar  coupling  could  arise  in  a  metal  with  the  nuclear  spins 
interacting  with  the  conduction  electrons.  The  existence 
of  such  a  coupling  can  best  be  explained  in  terms  of  a 
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simple  model. 

We  take  as  a  model  a  metallic  crystal  in  which  only 
two  of  the  nuclei  are  magnetic.  The  electron  wave  function 
in  the  periodic  lattice  will  be  scattered  by  the  interaction 
with  each  magnetic  nucleus .  The  effect  of  the  magne¬ 

tic  moment  of  a  nucleus  at  any  lattice  site  is  to  make  that 
site  favourable  for  an  electron  of  parallel  magnetic  moment  but 
unfavourable  for  an  electron  with  anti-parallel  moment.  In 

order  to  take  advantage  of  the  magnetic  interaction,  an 
electron  of  parallel  spin  will  distort  its  wave  function  so 
that  the  wave  function  is  slightly  larger  in  the  vicinity  of 
the  nucleus.  The  distortion  is  brought  about  by  mixing  in 
other  states  of  the  same  spin  orientation.  The  wavefunctions 
of  the  additional  Bloch  states  are  added  in  phase  with  the 
unperturbed  function  at  the  nucleus  in  order  to  interfere 
constructively  at  that  point,  but  because  of  the  spread  in 
wavelengths,  they  quickly  get  out  of  step  as  one  moves  away 
from  the  nucleus .  As  a  result  of  the  difference  between 
the  unperturbed  and  perturbed  functions ,  the  original  uni¬ 
form  distribution  of  spin  density  is  changed  to  have  a  damped 
oscillatory  behaviour.  The  total  electron  wavefunction  seen 
by  one  magnetic  nucleus  will  therefore  depend  upon  the 
spin  orientation  of  the  other,  thus  establishing  an  indirect 
spin-dependent  coupling  between  the  nuclei. 

The  strength  of  the  coupling  in  the  general  case 
may  be  calculated  from  the  total  electron-nuclear  Hamiltonian. 
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This  interaction  between  the  nuclei  and  the  electrons  £ 
is  written  as 


Hen  ?  ci  ii  *  *  -i  6  -iJ 


+  l  Di  li  *  l  hi  *  , 
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where 
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and 


(4) 


The  first  term  in  equation  (1)  is  the  'contact1  interaction 
between  the  electrons  and  nuclei,  whilst  the  second  term 
gives  the  dipolar  interaction.  Equation  (1)  may  be  rewritten 
in  the  form 


H 

en 


(5) 


where  G.  does  not  involve  the  nuclear  spin  coordinates. 

—l 

The  interaction  energy  associated  with  H  is  given,  to 
second  order,  by 


AE 
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The  terms  in  AE  such  as  <  0  |  lb  |  0  > 
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arise  if  each  nucleus  was  isolated  from  the  other  nuclei. 

This  energy  is  associated  with  parts  of  the  electron-nucleus 
interaction  which  give  rise  to  the  Knight  Shift  and  the 
relaxation  processes  in  metals.  The  cross  terms,  i.e. 
those  involving  two  nuclei  e.g.  <0|H^|n>  <n | H ^ | 0> (E^-E^  ^ , 

represent  the  interaction  energy  associated  with  the  electron- 
coupled  nuclear-nuclear  interaction.  These  may  be  expressed 
as 


(Oa|H.  |  not')  (na  '  |  H  .  1 0  a) 

AE  =  E  E  - - - ^ -  +  c.cv 

0a  i> j  na'  Ea+E  -E  -E  , 

J  0  a  n  a ' 


(7) 


where  |na')  is  a  product  of  a  many-electron  state  | n)  with 
energy  E^  and  the  many-nucleus  state  |  ot ' )  with  energy  E  ,  , 

|0a)  being  the  ground  state.  Substituting  equation  (5)  into 
equation  (7)  we  obtain 
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where  the  nuclear  energies  have  been  ignored  in  comparison 
with  the  electron  energies .  This  energy  AEq^  is  just  that 
which  would  be  expected  from  a  first-order  perturbation 
contribution  of  an  extra  term  in  the  nuclear  Hamiltonian 
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where  is  the  spin  coupling  tensor  between  the  i  and  j 

nuclei.  Taking  the  electron  wave  functions  to  be  products 
of  Bloch  functions,  it  is  possible  to  calculate  the 
(Abragam,  1961;  Slichter,  1963). 

The  most  important  term  arises  from  the  cross  terms 
involving  the  contact  interaction  between  the  nuclei  and 
the  electrons .  This  is  given  by 
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where 


(c.c.)  b.  .h  =  Z  |u.  ,  (0)  I  2|U.  (0)  I  2cos{  (k-k1)  .  (R.-R.)  } 
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U,  (0)  is  the  value  of  the  spatial  part  of  the  electron 
wave  function  at  the  nucleus  and  f (k)  is  the  Fermi  probability 
function.  It  can  be  seen  from  the  form  of  Hcontact_contact 
that  it  has  the  appearance  of  an  exchange  interaction. 
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Hence  its  name  of  'pseudo-exchange'  coupling. 

It  is  not  possible  to  evaluate  the  summation  with¬ 
out  either  some  further  approximations  or  some  explicit 
information  on  the  k  dependence  of  the  wave  functions  and 
energy.  We  assume  that  the  energy  surfaces  are  spherical 
and  there  is  an  effective  mass  m*(k) .  As  the  greatest 

contribution  to  J .  .  comes  when  k  -  k'  -  k^  .  lu,  (0)  I 

f  '  1  k  1 

2  — 
and  |u,  ,(0) |  may  be  replaced  by  a  value  appropriate  to  k 

and  k '  near  the  Fermi  surface  and  the  sums  may  then  be 

evaluated  (Slichter,  1963)  .  This  leads  to  a  value  of  the 

scalar  coupling  constant  Jh  ^  given  by 

h  Jij  =  "  If  Ye2  Yi  Yjldm*  |Ukf(0)  |4kf4F(2kf|R.j  |)  (12) 

where 

tp / \  _  sin  x  -  x  cos  x 

F(x)  -  - 7 -  .  (13) 

X 

This  is  the  Ruderman  -Kittel  approximation  (Ruderman  and 

Kittel,  1954) .  The  oscillatory  part  F(x)  is  plotted  in 

figure  (3.1)  .  The  dependence  on  |  U,  (0)  shows  that  the 

kf 

coupling  will  be  greatest  for  a  nucleus  with  a  large  Z. 

The  second  term  from  H  __  arises  from  the  cross - 

eft 

terms  between  the  contact  and  dipolar  interactions . 
is  given  by 

Hcontact-dipolar  ~  +  +  CC 
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Figure  (3.1) 

The  spatial  dependence  of  the  Ruderman-Kittel  'pseudo- 

exchange'  interaction  F  (x) .  The  figure  also  shows  the 

distances  of  the  neighbouring  nuclei  in  white  tin. 

(x  ) ^  indicates  the  value  of  x  taken  for  the  start 
o  Sn 

of  the  continuum,  as  discussed  in  chapter  4. 
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where 


B.  . 
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c .  D  .  E  E 
1  3  n  4 
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x  (n|St  .  .  IjlO)  (VV"1  (14) 

Summing  over  the  electron  spins  as  previously,  we  obtain  a 

traceless  coupling  (Abragam,  1961)  1^  .  J^j  .  1^  where 

JV.  is  proportional  to  P.  ,  +  P..  +  cc  with 
•  J  13  3 1 

E±i  =  E  (k|5(rk-R^|k')  <k'  |A  |k)  (Ek-Ek,)_1  .  (15) 

J  kk '  —  J  —  — 


This  coupling  will  be  removed  by  the  macroscopic  rotation 
(see  chapter  2) .  In  special  cases ,  (Bloembergen  and 
Rowland,  1955)  it  may  be  shown  that  this  coupling  takes 
the  form  of  a  pseudo-dipolar  interaction 
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but  there  is  no  group  theoretical  argument  that  this  is 
always  the  case  (Abragam,  1961)  . 

The  final  terms  in  involve  the  cross-products 

between  the  electron-nuclear  dipole  interaction  terms . 

This  may  be  expressed  as  a  sum  of  two  tensors, one  with  the 

II 

form  of  J.'.  1  and  the  other  with  the  form  J.  .  .  These 

13  =  =13 

terms  are  much  smaller  than  those  already  considered  and 
may  be  thought  of  as  correction  terms  (Abragam,  1961) . 
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CHAPTER  4 .  TIN 


A)  Theoretical  analysis 

Van  Vleck  (1948)  showed  that  spin-spin  coupling 
between  nuclei,  corresponding  to  a  Hamiltonian  of  the  form 


H 


ss 


J  h  I1  .  I  , 


has  two  effects  on  the  resonance  line-shape.  If  the  coup¬ 
ling  is  between  unlike  nuclei,  the  second  moment  is 
increased,  implying  that  the  line  is  broadened.  When  the 
nuclei  are  alike,  however,  the  coupling  produces  an  increase 
in  the  fourth  moment  but  leaves  the  second  moment  unchanged. 
This  means  the  line  narrows  but  there  is  an  intensity  shift 
towards  the  wings.  A  number  of  studies  have  been  made  on 
tin  spectra  using  the  ’method  of  moments '  (Karimov  and 
Schegolev,  1961;  McLachlan,  1968;  Schone,  1969;  Sharma 
et  al.,  1969;  Alloul  and  Deltour,  1969).  This  approach 
is  difficult.  The  experimental  line  shape  is  somewhat 
Lorentzian,  implying  that  there  is  considerable  intensity 
in  the  wings  which  makes  accurate  determination  of  moments 
difficult  (Sharma  et  al.,  1969) .  In  addition,  interpreta¬ 
tion  of  the  results  is  complicated  by  the  necessity  of 
unravelling  the  effects  of  the  spin-spin  scalar  coupling 
from  those  of  the  spin-spin  'pseudo-dipolar '  coupling. 

A  second  approach  which  has  been  used  to  determine 
the  scalar  coupling  is  to  apply  the  Anderson-Weiss  theory 


* 
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(Anderson  and  Weiss,  1953).  This,  being  useful  when  the 
line  is  very  exchange-narrowed,  is  only  applicable  to 
resonances  from  isotopically  enriched  tin.  The  experimental 
line  width  <5 to  is  given  in  this  theory  by 

6 co  =  tt<Ao)2>//3  J  (1) 

where  <Aw2>  is  the  theoretical  second  moment  in  the  absence 

of  exchange  narrowing  and  J  is  an  average  value  of  the 

scalar  coupling  constants  over  the  nearest  neighbours 

(McLachlan,  1968;  Sharma  et  al.,  1969). 

In  this  thesis  an  entirely  different  approach  is 

used.  An  attempt  is  made  to  build  up  the  tin  spectrum 

from  first  principles .  This  is  possible  in  natural  tin 

because  of  the  low  isotopic  abundance  of  the  magnetically 

n  7  1 1  q 

active  nuclei  (Sn  ,  7.67%;  Sn  ,  8.68%)  and  the  fact  that 
the  macroscopic  rotation  has  given  the  Hamiltonian  the 
same  form  as  for  a  liquid.  This  allows  the  possibility 
of  applying  high  resolution  nmr  spectroscopy  techniques 
in  a  situation  where  they  would  not  normally  prove  useful. 
The  low  isotopic  abundance  permits  the  splitting  of  the 
sample  into  groups  of  interacting  nuclei;  each  group  being 
effectively  isolated  from  the  next.  The  distances  from  a 
nucleus  to  its  four  shells  of  nearest  neighbours  are 
3.016  8,  3.175  A,  3.76C)8  and  4.1408  (see  figure  (4.1)). 

From  an  examination  of  the  form  of  the  Ruderman-  Kittel 
interaction  (figure  (3.1))  we  see  that  after  the  third 


Figure  (4.1) 


Diagram  of  the  tin  crystal  structure. 
3  and  4  indicate  the  nearest  neighbou 
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nearest  neighbour,  the  spin-spin  coupling  is  small.  We 
therefore  consider  in  some  detail  the  effects  on  a  central 
nucleus  of  its  neighbours  in  the  first  three  shells  and 
treat  the  other  nuclei  as  part  of  a  continuum.  The 
analysis  is  facilitated  by  the  fact  that  Sn117  and  Sn11^ 
both  have  spin  However,  provided  the  macroscopic 

rotation  has  removed  any  effect  of  quadrupolar  coupling, 
the  technique  could  easily  be  extended  to  cover  spectra 
involving  nuclei  with  I  >  ^-. 

The  nuclei  interacting  with  the  central  nucleus 
are  situated  in  shells  having  4,  2  and  4  sites  respectively. 
We  can  determine  the  probability  of  finding  certain  groupings 
of  nuclei  in  the  sample.  These  are  given  in  table  I.  With 
techniques  similar  to  High  Resolution  nmr  it  is  possible 
to  calculate  the  spectra  that  arise  from  these  configurations. 
In  metallic  tin,  the  presence  of  the  electron-nuclear 
interaction  very  effectively  relaxes  the  nuclei;  this 
broadens  the  resonance  and  leaves  little  possibility  of 
resolving  individual  lines. 

Taking  one  of  the  isolated  groups  we  construct  an 
ensemble  of  systems  each  with  the  Hamiltonian  H  of  the 
group , 

H  =  E  +  +  G  ( t)  +  H  -  ( t)  (2) 

o  1  rr 

The  time  independent  term  Eq  is  the  Hamiltonian  of  the 
isolated  nuclei.  It  has  eigen-functions  ^ ^  , 


etc . 
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Table  I 

Probability  of  finding  certain  configurations  of  tin  nuclei 
in  a  sample  of  natural  tin. 
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Probabilities  Are  Divided  Amongst  The  Various  Shells  As  Follows: 
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with  eigenvalues  a',  3'/  etc.  The  time-independent  term 
E^  is  produced  by  the  electron-nuclear  interactions  and 
is  responsible  for  the  Knight  Shift  and  the  electron-^ 
induced  spin-spin  coupling.  By  use  of  perturbation  theory 
it  is  possible  to  determine  the  eigen- functions  \p&,  ip^, 
etc.  of  the  Hamiltonian  Eq  +  E^.  The  Hermitian  perturbation 
G(t)  is  random  in  time.  It  takes  into  account  the  effect 
on  the  nuclei  of  the  random  magnetic  fields  produced  by 
the  electrons.  G(t)  is  therefore  responsible  for  the 
electron  part  of  the  nuclear  relaxation.  The  externally 
applied  rf  field  Hamiltonian  is  denoted  by  H^ft)  •  In 
general  the  wave  functions  of  the  ensemble  can  be  written 

Tp  =  E  C  \p  (3) 

—  ora  ■ 

a 

We  define  the  density  matrix  in  the  Schrodinger  represen¬ 
tation  by 


(4) 


where  <>  indicates  an  ensemble  average.  The  density  ma¬ 
trix  elements  satisfy  the  following  differential  equation 
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where  pftftl(T)  is  the  thermal  equilibrium  value  of  pQQI 


$3 

(Redfield,  1957).  The  R's  are  defined  by 


33 


2hR  t00i  ~  J  o  ini (a,"3l)  +  J  0  i0,(a-3) 
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The  J's  are  the  random- field  spectral  densities  and  are 
defined  by 


+  C0 


Ja3a'3'(w)  =  /  <  (a  |  G(t)  |  3)  (3  '  |  G  (t+x)  |  a  ')>exp  (-iwx)  dx  .  (7) 


It  may  be  shown  (Redfield,  1957)  thatR^i^i^  0  unless 


(E  -E  ,-E0+E0 ,)  t  <<  h,  (8) 

a  a '  3  3  c 

where  t  is  the  correlation  time  of  the  random  fields .  We 
c 

may  use  the  density  matrix  in  determining  the  expectation  value 
of  an  operator  M.  It  is  given  by 


<M>  = 


P  ,  M  ,  „ 

'act'  a'a  , 


(9) 


1  a , a ' , a" 

Having  introduced  the  density  matrix  we  make  use  of  it 
in  determing  the  tin  spectrum. 


Case  1:  No  active  nuclei  in  next  three  shells. 


In  this  instance  there  is  no  spin-spin  coupling 
present  and  the  spectrum  is  a  single  line  whose  width  is 
by  relaxation  broadening.  We  shall  however 
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treat  this  simple  case  in  some  detail  to  show  the  approach 
to  be  used  in  the  more  complicated  systems. 

The  Hamiltonian  of  the  single  nucleus  is 


H  =  -y'h  I  (1-K)  H  +  G  ( t )  +  H  .(t) 
i  o  rf 


(10) 


where  K  is  the  isotropic  Knight  Shift  constant.  The 
isolated  nucleus  has  two  energy  states,  with  its  spin 
either  parallel  to  the  magnetic  field  (state  | 1)  )  or 

antiparallel  (state  | 2) ) .  The  density  matrix  therefore 

has  four  elements  p^,  £>22'  pi2'  p21*  T^e  field 
produces  transitions  between  states  |  1)  and  |  2)  .  Writing 


(1|  Hrf(t)  I  2)  =  H12  , 
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we  have  from  equation  (5) 
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In  the  limit  of  small  rf  fields,  i.e.  ignoring  terms  of 
2 

the  order  we  ^ave  f rom  equation  (12) , 


11 

Pll(T) 

(15a) 

22 

p22 ^  • 

(15b) 

Taking  E 2  >  E^  and  assuming  that  H  ^(t)  is  sinusoidal  at  a 
frequency  oo,  we  try  a  solution  of  the  form 


Taking  H 
we  have 


iwt 

12 

r  e 

12 

-icot 

21 

r216 

rf(t)  =  V  cos  cot  =  V(exp(icot) 


+  exp (-imt) ) /2 


(16a) 

(16b) 

(17) 


r12lt0  <R1212  +  1Wo)rl2  +  V12(pll  p22>/2 


(18) 


where  E^  -  E_  -  "H  co  .  To  evaluate  the  R's  it  is  necessary 
2  1  o  J 

to  make  a  number  of  approximations.  First  is  that  the 
random  fields  produced  at  the  nuclear  sites  by  the  elec¬ 
trons  do  not  have  any  preferred  direction.  (This 
assumption  will  be  discussed  in  detail  later.)  We  also 
assume  that  the  frequency  spectrum  of  the  random  fields  is 
white  to  frequencies  much  greater  than  the  Larmor 
frequencies  of  the  nucleus .  This  implies  that  inequality 
(8)  is  satisfied.  As 


(E  -Ej  /h  -  10^  sec  "^and  t  -  10 
a  3  c 


sec 


. 


•  • 
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(t  being  the  electron  spin-lattice  relaxation 
time),  this  is  not  unreasonable.  We  can  now  evaluate 
the  R's.  From  equation  (6) 


2h2R 


1212 


=  2<J1122(0)  " 


JHll(°)) 


J2121(w) 


J1212 


J2222 


(19) 


Now  G  ( t)  =  -  v  "h  H  (t)  .  I 

'  n  r 


(20) 


where  is  the  random  magnetic  field  produced  by  the 
electrons  at  the  nuclear  site,  so  that 


+  00 

J  ,  ,  (to)  =  /<(ct|H  (t)  ,l|e)  (e1  |H  (t+T)  .I|a')  > 

o t  a  p  j-  j- 


2  2 

x  y  -ft  exp (-iwx) dx 


=  S  2h2Y  2  (a  1 1  I  6)  <  e  •  1 1  ,|a')k  ,  (to)  (21) 

,  n  lql  lq  l  qq 

where 


+  00 

2k  ,  (w)  =  /  <H  ( t)  H  (t+x)  >  exp  (-ioox)  dx  (22) 

-«  q  q ' 


and  H  (t)  is  the  qth  component  of  H  (t) .  If  the  random 

rq  — 

field  in  perpendicular  directions  are  assumed  to  be  un¬ 
correlated  then 


<  H  ( t)  H  (t+x)  >  =  0  (23) 

q  q' 


unless  q'  =  q.  We  have  that 


$  m  :I  SHU  im 
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=  -J1122(l0)  =  Yn2-ft2kz:>>/2 

and  J2121((o)  =  Yn?«2(kxx(“)  +  kyy<“>>/2  • 

The  assumption  of  a  white  frequency  spectrum  for  the 
random  fields  leads  to 


(24) 


(25) 


k  (to)  =  k  (0) 

qq  qq 


(26) 


If  the  random  fields  have  values  independent  of  direction 
then 


and 


k  =  k  -  k 
xx  yy  zz 


R1212 


+  k  ) 

yy 


(27) 


=  -2W 


(28) 


where  W  =  Yn^k  .  Substituting  equation  (28)  back  into 
equation  (18)  we  obtain 

r12  =  “^-Vi2  ^ pn” P12 ^  (wo~w)  -  2W } 2h ]  1.  (29) 


Now 


P1J  =  exP  (“Ej^/kT)  {exp  (-E^/kT)  +  exp(-E2/kT)} 


-1 


-  {exp(-E1/kT)  }/2  (30) 

so  that  pn-p22  ~  ~fr^0/2kT  .  (31) 

Using  the  density  matrix  we  find  that  the  magnetization  in 
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the  x  direction  is 

,  / 

<Mx(t)>  =  P12  (2  |  Mx|  1)  +  p21(l'|Mx|2) 

=  r12  (2  |  Mx|  1)  exp  (icot)  +  r21(l|.M  |  2)  exp(-iwt) 

=  2  Real  [r12  exp  (icot)  ( 2  |  Mj  1)  ]f  (32a) 

which  by  definition  equals  the  real  part  of 

X(w)  Hxoexp(iwt).  (32b) 

Therefore,  we  have 

X<<*>)  =  2r12(2|  Mxl  D/Hxo,  (33) 

Using  the  fact  that 


V,  0  =  -  ( l|  M  H  |  2) 
12  1  x  xo1 


(34) 


leads  to  a  value  of  the  nuclear  susceptibility 


X(to)  =  i  |  (1  |M  |  2)  |  2to  {2kT  (i  { oj  —  oj)  -2W)  }  1 

X  u  u 


=  4i  Z  { i  ( co  —  <o)  -  2W} 
o  o 


-1 


(35) 


2  2  - 1 

where  Z  is  defined  as  Y  fi  co  (32kT)  .  We  see  that  the 

o  no 

line  is ,  as  expected,  a  Lorentzian  with  a  width  proportional 


to  W. 


. 


Case  2: 


One  active  nucleus  in  next  three  shells . 
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We  now  suppose  that  there  is  an  other  magnetically- 
active  nucleus  within  the  first  three  shells  .  In  this 
case  we  would  expect  two  types  of  spectra;  one  when  the 
two  nuclei  are  alike,  the  other  when  they  are  not.  In 
the  first  instance,  the  presence  of  the  spin-spin  coupling 
does  not  change  the  spectrum  from  that  of  an  isolated 
nucleus,  but  in  the  other  case,  the  system  produces  what 
a  chemist  would  call  an  AX  spectrum.  As  the  nuclei  are 
unlike^  their  resonance  lines  do  not  overlap  and  can  be 
treated  separately.  The  effect  of  the  J  coupling  of  the 
A  nucleus  on  the  X  spectrum  is  to  split  it  into  two  lines . 
It  is  as  though  the  X  nucleus  would  be  found  in  two  magne¬ 
tic  environments;  one  when  the  A  nucleus  has  spin  parallel 
to  the  main  field,  and  the  other  when  antiparallel.  The  J 
coupling  from  A  acts  as  an  'extra'  magnetic  field  at  the 
X  site.  The  random  fields  at  the  A  nuclear  site  can  be 
thought  of  as  flipping  it  from  spin  up  to  spin  down  and 
vice-versa.  This  change  effectively  switches  the  magnetic 
environment  of  nucleus  X.  If  the  flip  rate  is  high 
enough ,  the  effect  is  to  decouple  the  X  nucleus  from  the 
A  so  that  the  X  doublet  collapses  into  a  line  centered  at 
the  natural  frequency  of  X.  The  coalescing  of  the  lines 
is  a  common  enough  effect  in  high  resolution  nmr 
(Arnold,  19  56)  .  There, it  is  often  produced  by  the  jumping 
of  nucleus  X  from  a  molecule  where  the  A  spin  is  up  to 
another  where  it  is  down.  This  also  has  the  effect  of 
changing  the  X  magnetic  environment.  In  liquids  the 


- 

- 


WS  <v  fif  3 
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collapse  of  the  lines  is  very  apparent,  but  in  a  metal 
the  effect  of  the  relaxation  broadening  (determined  by  W) 
is  so  great  that  the  multiplets  will  not  be  resolved  and 
their  collapse  will  be  manifested  only  by  a  slight 
narrowing  of  the  line. 

It  should  be  realized  that  there  is  a  difference 
between  the  disappearance  of  the  satellites  under  this 
flipping  process  and  the  disappearance  of  the  dipolar 
broadening  when  the  sample  is  spun.  Rotation  keeps  the 
relative  spin  positions  unchanged  and  produces  a  spatial 
average  leaving  the  multiplet  structure.  The  flipping 
however  changes  the  spin  positions  and  thus  collapses 
the  multiplet. 

We  shall  now  proceed  to  show  that  when  the  two 
nuclei  are  identical  there  is  only  one  line.  The 
Hamiltonian  has  the  form 

H  =  -y  *hl  (1-KjH  -  Y^.111  (1-K  ,)H 

A  z,  A  o  A'  z_  ,  A '  o 
A  A ' 

+  JAA'hIA  *  1 a '  +  G(t)  +  Hrf (t)  *  (36) 

The  time-independent  part  of  the  Hamiltonian  has  four 
eigenfunctions.  Three  of  these  are  associated  with  a  total 
spin  of  1,  and  the  other  is  of  spin  0.  It  is  not  possible 
for  the  applied  rf  field  to  induce  transitions  between 
states  of  different  total  spin.  The  nmr  lines  will  there¬ 
fore  only  be  produced  by  transitions  between  the  states 
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of  total  spin  1  corresponding  to  m  =  -1,  0  and  +1  where 
mz  is  the  z  component  of  the  total  spin.  It  can  be  shown 
that  these  states  have  energies  y  H  ti  +  J  ,h/4, 

J  ,h/4,  and  -y  H  *K  +  J  ,h/4  respectively.  It  follows 

O  AA 

that  the  energy  required  to  induce  transitions  is  indepen¬ 
dent  of  J  ,  and  the  spectrum  is  the  same  as  for  the 
isolated  nucleus. 

For  non-identical  nuclei  the  Hamiltonian  becomes 

H  =  -Y*hl  (1-KJH  -  Y^h  I  (1  -  K  )H 
'A  z  A  o  X  z  X  o 

X 

+  J^h  IA  .  Ix  +  G  ( t)  +  Hrf  (t)  .  (37) 

The  eigenfunctions  of  the  time-independent  part  of  the 
Hamiltonian  are  |++) ,  |+-) ,  |-+)  and  |--)  which  we 

number  |l),  \2)t  |3)  and  | 4)  respectively.  The  energy 

level  scheme  is  shown  in  figure  (4.2)  where  we  have  chosen 
|e4  -  E2 |  >  |e4  -  E3 | ,  without  loss  of  generality.  The 

presence  of  the  random  fields  produces  transitions  between 
the  eigenfunctions.  We  take  the  simplest  form  of  G(t) 
which  is  a  sum  of  'single-nucleus '  operators 

G(t)  =  -ti(YAHA(t)  •  I  +  YxHx(t)  .  Ix)  (38) 

where  H,  and  H„  are  the  random  fields  at  the  A  and  X 
A  X 

nuclear  sites.  We  shall  assume  that  the  random  fields  at 
the  two  sites  are  independent.  From  the  form  of  G(t)  we 


f. 


. 
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Figure  (4.2) 

Energy  level  diagram  for  the  AX  type  grouping  of  nuclei 
showing  the  allowed  transitions. 
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note  that  will  induce  transitions  between  |4)  and  | 2) 
and-  also  |  3)  and  |  1)  whilst  H  induces  transitions  between 

i  ^ 

| 4)  and  | 3)  and  also  |2)  and  | 1) .  Imagine  that  we  are 
looking  at  the  X  spectrum^  then  the  quantity  of  interest  is 
the  x  component  of  the  magnetization  of  the  X  nucleus . 

This  is 


(t)>X  =  Pl2(2'Mx|l)  + 
+  p  3  4 ( 4 | |  3 }  + 

=  2  Real{p12  (2|  Mx|  1) 


P21U|Mx|2) 

p43(3|mxU) 

+  P 3 4  <  4 1  |  3)  }  . 


(38) 


Terms  such  as  (1|M  |3)  deal  with  lines  of  the  A  spectrum 
and  will  be  zero  near  the  X  lines .  Making  use  of  equations 
(5) ,  (6)  and  (7)  we  obtain,  in  the  low  rf  field  approxi¬ 

mation, 


pii = 

pii<T>  -  1 

=  1  +  4, 

(39) 

dp12 

dt 

lw12p12  + 

^P11_P22^  H21^ 

+  R1212p12 

+  R1234p34  ' 

(40) 

dp34 

dt 

=  lw34p34  + 

^p33_p44^H43//h 

+  R3  4  3  4  p  3  4 

+  R3412p12 

(41) 

=  to  -  Jaxtt  and  w34  =  with  wx  being 


where  w 


' 
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the  natural  frequency  of  the  X  nucleus . 
R's  leads  to  the  equations 


dp12 

“dt  lt012P12  +  l(°ll  p22)H21/,h 


-  P12(2Wx  +  WA)  +  P34W 


and 

dp34 

dt 


lw34p34  +  l(p33  p44)H43//h 


"  p34(2WX  +  V  + 


P12WA' 


Evaluating  the 


(42) 


(43) 


The  W's  are  as  previously  defined  with  the  suffixes  denoting 

the  nuclei  with  which  the  random  fields  are  associated  with. 

Ignoring  terms  involving  W^,  the  equations  for  p^^  an(^  P34 

look  just  like  those  for  the  density  matrix  of  a  two-spin 

system  with  the  nuclei  isolated  (see  equation  (13)).  The 

form  of  the  lines  is  Lorentzian  with  a  linewidth  determined 

by  2W  .  The  coupling  terms,  which  involve  W  ,  produce  the 
X.  A 

collapsing  of  the  lines,  a  phenomenon  which  might  well  be 
called  'relaxation  narrowing'.  Trying  a  solution  similar 
to  that  for  the  single  nucleus  we  take 


H^^(t)  =  V  cos  cot 


p12 

P34 


12 


34 


exp ( iwt) 
exp ( iwt) 


which  leads  to 


. 


. 
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X12r12  +  W  r34  lV12(pll  p22^2h 


and 

W  r12 

+  X34r34  lV34^p33  p44^ //2h 

(44b) 

where 

X.  .  = 
1 3 

i  ( w .  .  -co)  -  3W  . 

13 

(44c) 

Now 

p .  .  = 

'll 

4 

exp(-E./kT)  /  l  exp(-E./kT) 

j=i  3 

(45a) 

so  that 

p .  .  - 

'll 

pi+l,i+l  “  fi“o/4kT- 

(45b) 

From  equation 

(44)  we  obtain 

ii 

CO 

u 

2iZoHxo(X12-“)(X12X34-w2)_1(Y«)"1 

(46a) 

and 

r12 

2iZoHxo(X34-“>  (X12>.34-W2)-1(Y«)-1. 

(46b) 

We  therefore  have 

=  2(r12(2|Mx|l)  +  r34(4|Mx|3))Hxo-1 

=  2iZo(A34  +  X12  -  2W) (X12X34  -  W2)-1  .  (47) 

If  W  is  small  compared  to  2tfj>  equation  (47)  describes  a 
spectrum  consisting  of  lines  at  co  +  uj  broadened  by  W. 

If  2ttJ  is  small  compared  to  W,  equation  (47)  describes  a 
spectrum  with  a  line  at  oo  broadened  by  J  and  W. 

A 


f 


a 


4 

r 


. 
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Case  3:  Two  active  nuclei  in  next  three  shells 


When  both  nuclei  are  alike  and  identical  to  the 


central  nucleus ,  we  again  find  that  the  nmr  spectrum  of 

the  group  is  a  single  line.  The  spectrum  is^  however^  split 

into  a  total  of  15  lines  when  only  two  of  the  nuclei  are 

the  same  isotope.  It  is  possible  to  gain  an  insight  on 

how  this  spectrum  arises  by  considering  the  case  where  all 

the  nuclei  are  very  different.  In  this,  so  called  AMX  spectrum, 

the  line  from  the  A  nucleus  is  split  into  two  by  the  J 


coupling,  and  each  of  these  lines  is  further  split  by  the 

coupling.  A  similar  thing  happens  to  the  M  and  X  lines, 
giving  rise  to  a  total  of  12  lines  .  As  we  let  the  A  and 
M  nuclei  become  more  alike>  the  A  and  M  spectrum  come 


AA'  part  of  the  AA'X  spectrum.  The  four  lines  near  w 

A 

become  six  by  the  appearance  of  two  other  lines  associated 
with  the  simultaneous  flipping  of  the  A  and  A'  nuclei. 

The  fifteenth  line  is  a  combination  line  (near  C0A+wAi-a)x^ 
involving  A,  A'  and  X  nuclei  and  has  zero  intensity  unless 
X  is  similar  to  A. 


The  Hamiltonian  for  the  three  nuclei  A,  A'  and  X  is 


H  -  -h  E  y .  I 
i=A , A ' ,X 


.  I  ( 1-K  . )  H  +  G(t) 
l  z .  1  o 

l 


(48) 
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The  eigen-functions  of  the  time-independent  part  of  the 
Hamiltonian  are  given  in  table  II.  The  determination  of 
the  eigen-functions  is  greatly  simplified  by  the  fact  that 
the  basis  product  functions  with  different  X  spins  cannot 
be  mixed  .  Figure  (4.3)  shows  schematically  the  various 
allowed  transitions  produced  by  an  external  rf  field.  The 
transition  frequencies  have  been  given  by  Pople,  Schneider 
and  Bernstein  (1959)  and  are  shown  in  table  III. 

Let  us  first  imagine  that  we  are  considering  the 
X  part  of  the  spectrum.  This  means  that  the  X  atom  is  at 
the  centre,  surrounded  by  the  ten  sites.  The  lattice 
structure  in  tin  is  such  that  the  probability  of  finding 
the  A'  and  A  nuclei  closer  together  than  the  fourth-neares t 
neighbour  is  small.  We  can  therefore  simplify  the  calcu¬ 
lation  by  assuming  that  J^,  is  zero.  This  means  the  X 
spectrum  now  only  consists  of  four  lines.  In  this  case  we 
have 

HxoX<“>  =  2  ?.eij(jlMxli)  (49) 

1J  J 

where  ij  take  the  values  12,  36,  45  and  78.  Proceeding 
as  before,  and  taking 


G(t)  =  -  E  y,  "h  H  (t)  .1  , 

k=A , A ’ ,X  —  — 


(50) 


we  find  that  r^/  r36'  r45  and  r78  satisfy  the  simultaneous 
equation  (with  coefficients  given  in  matrix  form) 


■ 


Table  II 


Wave  functions  for  the  AA*X  type  grouping  of  nuclei. 


Wave  Functions  For  AAX  Spectrum 


1 1  )  =  I  +  +  +  ) 


17  )=!--+} 


I  2  )=  I+  +  -) 


|8)=  I - ) 


|3)=cos0  I+-  -  +)  4-sinO  I-  +  +  ) 


4 


I  4  )--sin()  Ih —  4- )  A  cos  0 1 — FT) 
T+  + 


|5  )=  cos  0  1  +  — )  4-  sinO  I-  +  -) 


|  6  )  =  -sinO  l+-  -  )  4  cos  0  I-  +  ~) 

T-  T  + 


whe  re 


D  cos2(p+  =  (J2-  J3)/4 


D+sin  2(J)±  =  J / 2 


2  2  1/2 

D±=((J2-J3)  /4+J,  )  / 2 
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Figure  (4.3) 

Energy  level  diagram  for  the  AA'X  grouping  of  nuclei  showing 
the  allowed  and  forbidden  transitions. 


STATE 


X  LINES 
AA'  LINES 

FORBIDDEN  COMBINATION  LINE 


62 


Table  III 


Transition  frequencies  for  the  AA'X  grouping  of  nuclei. 


Transition 


FREQUENCIES  FOR  AA'  SPECTRUM 


^68 

CUA 

^47 

WA 

cu 

25 

TZ 

loa 

^13 

— 

WA 

^58 

— 

^A 

^37 

WA 

^26 

WA 

^14 

WA 

■ran  SIT lOi 

UJ 

78 

— 

UJ 

X 

^35 

— 

UJ.n 

46 

— 

wx 

UJ^2 

— 

wx 

W45 

O) 

X 

W31 

— 

wx 

-  (2J.,  +  J„  +  J3)/4-  -  D_ 

-  (2 J  -J  +J  )/4  -  D 

12  3  -t- 

+  (2 J  “J  +  J  )/4  -  D 

12  3 

+  (24,  +J2+J3)/4  - 

-  (2  J  +J  +J  )/4  +  D 

12  3 

-  ( 2  -  J^+ )/4  +  D  + 

+  ( 2 J  -J  -J  )/ 4  +  D 

12  3 

4  (2 J  +J  +J  )/4  4  D 

12  3  f 


-  (J  4  J  )/ 2 

2  3 

+  D  -  D 

4 

-  D  +  D 

4  _ 

t(J2  +  J)/2 

-D-D 

+• 


4  D+  4  D 


Forbidden  Line 


^27 


Spectrum 


- 
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y 


12 

W 


y 


-w 

o 


w 

36 

0 

W 


where  y .  .  = 
ij 

that 


iZo 

Yh“ 


(51) 


Using  equation  (51)  we  have 


X(o3)  =  -2iZQ(i  [to  -w]  -  4W) 

X  [4W{i(cox-co)  -  4W}  -  y45u36  -  P12y78] 

x  (u12y78V45y36  “  4W2{i(ux-lo)  -  4W})71  (52) 

We  now  examine  the  AA'  part  of  the  AA'X  spectrum. 
This  is  in  general  asymmetric  about  w  .  If,  however,  we 
imagine  that  the  A  nucleus  is  at  the  centre,  then  A'  and 
X  are  so  far  apart  that  is  zero.  This  has  the  effect 

of  symmetrizing  the  spectrum.  Proceeding  as  before  we  have 
that 


X  (w)  H 


xo 


where  nq  take  the 
Setting 


values  68,  47, 


25, 


(53) 

13,  58,  37,  26  and  14. 


r  exp(iojt) 
nq  * 


(54) 


we  find  that  the  r  's  satisfy  the  simultaneous  equations 

nq 

shown  in  table  IV  (case  1)  .  There  is  nothing  to  be  gained 


*' 
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Table  ET 
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LU 

X 

X 


LU 

Z 

I 

x 

lu 

X 

XI 
Q 


X 

O 

X 

$ 

00 

z 
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x 
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by  attempting  to  solve  this  equation  algebraically  and  the 
solution  is  therefore  obtained  by  computation  at  each  point 
in  the  spectrum.  The  problem  is  not  difficult  as  the  matrix 
of  the  coefficients  of  the  r's  has  a  number  of  special 
properties.  The  method  of  analysis  is  given  in  appendix  I. 

Case  4:  The  effect  of  more  distant  neighbours. 

In  the  previous  analyses  of  tin  data  (McLachlan,  1968; 
Alloul  and  Deltour  1969  ;  Sharma  et  al. ,  1969)  the  effect 
of  far  distant  neighbours  on  the  spectrum  has  been  ignored. 
Approaching  the  analysis  from  the  'high  resolution  nmr ' 
side  shows  that  the  distant  neighbours  may  actually  have  a 
pronounced  effect  upon  the  line  shape. 

We  may  approximate  the  effect  of  an  X  nucleus  in 
the  'continuum'  on  an  A  nucleus  at  the  centre  of  a  group  by 
examining  the  AX  spectrum  in  the  limit  of  high  flip  rates 
(i.e.  W  >>  2ttJ  )  .  We  find,  after  considerable  mathemati- 

nA 

cal  manipulation,  that* 

XA(co)  =  i4Zo(i{ioA-u}  -  2W')"1 
o  o 

where  2W'  =  2W  +  4tt  J  /2W.  The  2W  in  the  denominator  of 

1\2\ 

the  term  involving  J  arises  from  the  fact  that  the  X 
nucleus  is  changing  its  spin  direction  at  a  rate  2W.  Equation 

*1  would  like  to  thank  Dr.  McClung  of  the  Department  of 
Chemistry  for  bringing  this  expression  to  my  attention. 
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(57)  implies  that  the  multiplet  has  collapsed  into  a 
single  Lorentz  line  whose  width  is  slightly  greater  than 
the  normal  linewidth  because  of  the  presence  of  the  residual 
linewidth .  Each  interaction  with  an  unlike  nucleus  in  the 
continuum  increases  the  apparent  linewidth.  Summing  over 
the  continuum  we  have  that 


( 2W ' ) 


apparent 


=  2W  +  B 


2  2 

where  B  =  E4tt  J,,,/2W 

N  M 


(58) 

(59) 


with  the  summation  being  over  all  the  X  type  lattice  sites 
in  the  continuum.  Now 


E  =  /  2Tr2pr2J2  (r)  dr  (60) 

N  r 

o 

where  p  is  the  number  of  unlike  active  nuclei  per  unit 

volume  of  the  continuum  and  r  is  the  distance  to  the 

o 

start  of  the  continuum  from  the  A  lattice  site  (approximately 
4th  nearest  neighbour  distance) .  From  chapter  3  we  have 


2ttJ  (r)  =  A  F  ( 2kf r)  (61) 

so  that 

E  J^N  =  A2[X2TTkf3xo3)_:L  (62) 

N 

where  x  =  2k^r  .  Now  Jn ,  the  nearest  neighbour  spin-spin 
o  t  o  ± 

coupling  constant,  is  given  by 


2ttJ1  =  A  F(2kfrJ[) 


(63) 
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so  that 

^JAN  =  pJl  I12llkf  F2(2kfr1)x^]“^  (64) 

Taking  typical  values  for  the  constants  for  tin 

(kf  =  1.65  x  108cm  ”1,  F(2kfrx)  =  8.4  x  10_4,p  =  3.2  x  lo21 

-3  3 

nuclei  cm  ,  ^  =  150  ysec,  J1  =  3.1  x  10  Hz  and  xq  =  15). 

leads  to 


( 2W*  ) 


apparent 


2W  +  1050  sec 


(65) 


As  2W  is  approximately  6700  sec  ^  we  see  that  the  addition¬ 
al  width  from  the  far  neighbours  is  important  and  should  not 
be  ignored. 

We  shall  now  try  to  compute  the  effects  of  the 
like  nuclei  in  the  continuum.  Returning  for  a  moment  to  the 
simple  AA'  spectrum,  we  find  that  the  effect  of  J  ,  is  to 
make  the  A  and  A'  nuclei  exchange  their  spins  at  a  rate 
4ttJ  ,  .  If  the  A  nucleus  is  in  a  group  and  the  A'  nucleus 
is  in  the  continuum,  then  half  of  the  time  the  exchange  of 
spins  between  A  and  A'  leaves  A  in  a  different  magnetic 
environment  in  the  group  (i.e.  the  A  spin  flips) .  The  rest 
of  the  time  A  and  A*  exchange  equivalent  spins  so  that  the 
magnetic  environment  of  A  is  unchanged.  We  may  approximate 
the  effect  of  the  A'  on  the  A  as  an  increased  flip  rate  of 
A  .  The  actual  effect  of  the  A'  on  an  AX  spectrum  is  shown 
in  figure  (4.4)  .  The  upper  spectrum  shows  the  AA'  part 


67 


Figure  (4.4) 

Spectra  showing  the  effect  of  exchange  coupling  between  like 

nuclei . 

(a)  The  spectrum  for  a  group  of  A, A1  and  X  nuclei  when  there 
is  a  coupling  between  the  A  and  X  nuclei,  and  the  A' 
nucleus  is  isolated  from  the  others. 

(b)  The  same  spectrum  except  that  there  is  now  an  exchange 
coupling  between  the  A  and  A'  nuclei. 
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when  the  A'  is  not  coupled  to  the  A  nucleus.  The  A'  line 
is  not  shifted.  The  coupling  splits  the  A  line  into 

two.  Introducing  the  J  ,  coupling  narrows  the  spectrum. 

The  narrowing  could  be  brought  about  by  an  increased  jump 
rate  of  the  Aspin,  but  in  this  case  the  lines  would  not 
remain  sharp  but  would  become  broadened.  As  this  is  only 
an  approximation  of  what  is  actually  happening^  there  is , 
unfortunately,  no  theoretical  expression  equivalent  to 
equation  (57).  The  apparent  increased  flip  rate  must  there¬ 
fore  be  estimated.  Now  a  J  coupling  splits  the  A  line, 
of  half  width  2W,  by  an  amount  2ttJ  but  after  collapse 

the  spectrum  becomes,  from  equation  (57) ,  a  Lorentz  line 

2  2 

with  half-width  2W  +  4tt  J  /W.  We  can  approximate  this  line 
by  two  Lorentz  lines  of  linewidth  2W  separated  by  the  extra 
half-width,  i.e.  a  splitting  of  2tt2J^/W.  This  effectively 
is  the  same  as  cutting  the  strength  of  the  J  coupling  by  a 
fraction  ttJ/W.  The  effect  of  the  J  ,  coupling  was  to 
exchange  the  spins  at  a  rate  but  if  we  introduce 

the  effect  of  the  conduction  electrons  then  the  exchange 
rate  drops  to  4tt  J  ,/W.  We  can  now  use  this  expression  to 
determine  the  rate  AW  at  which  exchange  between  the  A  nucleus 
at  the  centre  and  the  A'  nuclear  spins  in  the  continuum 
occurs.  It  is  given  bv 

AW  =  l  2.72tt2/W,  (66) 

N 

where  the  factor  of  —  arises  as  only  half  of  the  A'  nuclei 


' 
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have  spins  such  that  their  exchange  with  A  changes  the 
A  environment.  It  follows  that 

AW  =  ti3PJ12[  6Wkf3F2(2kfr1)xo3]':}-  (67) 

If,  for  example,  we  include  the  effects  of  nuclei  in  the 
continuum  on  the  AX  spectrum  then  we  find  that 

(X12-M+B)r12  +  (W  +  AW)r34  =  2iZo/yh  (68a) 

(W  +  AW)r12  +  U34-/5-7+B)r34  =  2iZQ/yh.  (68b) 

The  extra  broadening  term  B  produces  an  increased  width, 
whilst  the  increased  jump  rate  AW  increases  the  terms  which 
collapse  the  multiplets  in  the  spectrum  (but  are  not 
included  in  the  terms  which  broaden  the  spectrum  ) . 

The  treatment  has,  however,  been  slightly 
inconsistent  as  the  nuclei  in  the  continuum  surrounding  A 
are  affected  i.e.  flipped  by  'their'  continuum  and  local 
neighbours.  We  must  therefore  make  the  expressions  for  B 
and  AW  self-consistent.  This  implies  that 

B'  =  B  2W(2W  +  AW')-1  (70a) 

AW'  =  ( AW) 2W ( 2W  +  AW')"1  +  AW"  (70b) 

where  AW"  is  the  effect/  on  the  nuclei  in  the  continuum, 
of  their  local  neighbours.  We  can  estimate  this  bv  taking 


. 


i 
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AW"  =  <Z  tt2J2(2W  +  AW')  1>t  ’ 

average  over 


(71) 


sum  over 
local  neighbours 


sample 


A  plausible  value  for  AW"  is  0.13  J^2(2W  +  AW')  ^ . 
Substituting  this  value  back  into  equation  (70b)  leads  to 


AW'  -  AW ( 1  +  0.07tt2J12W  2) 


(72a) 


and 


B  '  =  B  2W  ( 2W  +  AW  ( 1  +  0.0  7ti2J12W  2)  )  1 


(72b) 


We  can  now  use  these  expressions  to  correct  the  values  of 
X(oo)  for  the  effect  of  far  distant  neighbours  (cf  McConnell, 
1957) .  Defining 


W  =  W  +  AW' 

e  ? 


(73) 


then  the  isolated  nucleus  spectrum  becomes 


X(w)  =  i4Zo[i(a)o-a))  -  2W  -  B'l  1. 
The  X  part  of  the  AX  spectrum  becomes 


X.x(u)  =  2iZ„(*;4  -2W^)(X1,9  x;4  -W^  ) 


2V  -1 


o  34  12 


12  34  e 


(74) 


(75) 


where 

XV.  =  i  (^±  j“w)  "  2W  -  We  -  B'# 

Equation  (52)  describing  the  X  part  of  the  AA'X  spectrum 
must  be  modified  to  read 


X(«)  =  -2izoY(4weY  -  y45  y36 


P12  ^78  ' 


x  *y12  U78  y45  P36 


-  4W  2Y2)  1 

e 


(76) 


. 


. 
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where 

Y  =  i  (m  -to)  -  2W  -  2W  -  B' 
x  e 

(77) 

and 

y!  .  =  i(w.  .-to)  -  2W  -  2W  -  B' 

19  19  e 

(78) 

The  new  equations  for  the  AA*  part  of  the  AA'X  spectrum  have 
been  given  in  Table  IV  (case  2) . 

Case  5:  Other  configurations. 

These  account  for  the  remaining  20%  of  the  assemblies. 
If  the  nuclei  in  the  grouping  are  identical,  only  a  single 
line  is  expected.  In  other  cases  the  spectra  could  be 
calculated  but  the  proliferation  of  lines  to  be  considered 
makes  the  calculation  impractical.  Furthermore,  the  pro¬ 
bability  of  occurrence  of  each  configuration  is  small.  To 
ignore  them  entirely  would  be  incorrect,  so  instead,  we 
group  them  with  the  configurations  which  they  most  closely 
resemble.  Any  configuration  with  one  unlike  and  two  or  more 
like  nuclei  is  treated  as  if  it  consisted  of  one  like  and 
one  unlike  nucleus.  All  others  are  included  in  the  two  un¬ 
like  nuclei  group. 

Determination  of  a  value  for  W 

In  metals  at  room  temperature,  the  spin-lattice  re¬ 
laxation  time  is  mainly  determined  by  the  conduction 
electrons.  We  have,  for  an  isolated  nucleus  being  relaxed 
by  the  electrons  (Slichter,  1963)  ,  that 

T.  =  1/  (k  +  k  )  (79) 

1  '  v  zz  yy 

so  that  2W  is  equal  to  the  reciprocal  of  the  spin-lattice 


•  ■ 


r 
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relaxation  time.  The  assumption,  made  in  the  first  part 
of  this  chapter,  that 


k  =  k  =  k 
xx  yy  zz 


(80  ) 


implies  that  is  isotropic  and  this  can  be  tested  by 
measuring  the  orientation  dependence  of  T^.  This  has  been 
done  by  McLachlan  (1968)  who  finds  that  there  is  little 
variation.  The  absolute  value  of  his  T^  data  must  be 
considered  suspect  as  it  appears  that  when  compared 
to  other  work  there  is  evidence  of  some  systematic  error 
(Spokas  and  Slichter,  1959;  Dickson,  1969;  App,  1970). 
However,  a  systematic  error  should  not  affect  the  observa¬ 
tion  of  an  orientation  dependence  of  T^. 

There  are  a  number  of  values  of  T^  in  the  literature 
McLachlan  (1968)  together  with  Asayama  and  Itoh  (1962)  find 
values  of  T^  that  at  room  temperature  are  equivalent  to 


112  Msec. 


This  is  close  to  the  theoretical  value  of  T 


as  determined  using  the  known  tin  Knight  Shift  and  the 
Korringa  relationship  (Korringa,  1950).  However,  Spokas 
and  Slichter  (1959)  ,  Dickson  (1969)  /  Alloul  and  Deltour 
(1969)  and  App* (1970)  find  values  close  to  150  Msec. 

App's  work  was  actually  performed  using  the  same  crystal 


*1  would  like  to  thank  Dr.  Williams  of  the  Department  of 
Physics  at  U.B.C.  for  bringing  my  attention  to  this  work. 


73 


used  bv  McLachlan.  Dickson  (1969)  tested  two  other  samples 
of  tin,  one  of  high  purity  and  the  other  of  low  purity, 
and  obtained  T  's  consistent,  within  experimental  error, 
with  his  original  result.  We  shall,  therefore,  use  his 
value  of  in  the  determination  of  W. 


' 
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B)  Experimental  Results 

A  typical  absorption  derivative  nmr  signal  from  a  tin 
metal  sample  spun  at  5  kHz  at  the  'magic'  angle  is  shown  in 
figure  (4.5a).  Figure  (4.5b)  shows  the  resonance  of  an 
unspun  sample  to  the  same  scale.  The  complete  resonance 
of  the  unspun  sample,  shown  in  figure  (4.5c),  clearly  shows 
the  lineshape  asymmetry  produced  by  the  anisotropic  Knight 
Shift. 

The  lineshape  for  both  the  tin  isotopes  were  Lorentz- 

like  with  a  peak-to-peak  width  of  1.70  +0.05  kHz  for 

117  119 

Sn  and  1.90  +0.05  kHz  for  Sn  .  If  the  lineshapes  had 

been  totally  determined  by  the  spin-lattice  relaxation  times 

of  the  nuclei,  the  widths,  corresponding  to  T^'s  of  164 

ysec  and  150  ysec,  would  have  been  1.14  kHz  and  1.25  kHz 

117 

respectively.  The  Sn  nuclear  relaxation  time  was  cal¬ 
culated  using  the  fact  that  (T^) Sn117/ (T^) Sn119  =  Yggg/Yggy* 

It  was  found  that  there  was  no  change  in  the  resonance 
lineshape  on  increasing  the  rotation  speed  from  4.2  kHz  to 
5.7  kHz.  This  was  taken  as  evidence  that  the  rotational- 
satellite  lines  were  sufficiently  weak  and  far  from  the 
centre  line,  so  that  the  spectrum  could  be  described  by  the 

'liquid-like'  Hamiltonian  used  in  the  first  part  of  this 

\ 

chapter  (see  also  appendix  II) .  The  variation  of  the  reso¬ 
nance  linewidth,  as  the  turbine  was  tilted  relative  to  the 
magnetic  field,  gave  an  insensitive  determination  of  the 
'magic'  angle.  Instead,  the  disappearance  of  the  lineshape 
asymmetry,  associated  with  the  anisotropic  Knight  Shift/was  used. 


-- 
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Figure  (4.5) 

Experimental  Sn^^  resonances: 

(a)  Resonance  from  a  powder  sample  spun  at  5  kHz.  The 
resonance  frequency  is  6  MHz.  The  frequency  scale  is 

574  Hz/div.  The  resonance  was  obtained  by  making  2048  runs 
over  a  period  of  48  hrs . 

(b)  Resonance  from  an  unspun  sample  to  the  same  scale  as  (a) . 

(c)  A  resonance  from  an  unspun  sample  on  a  smaller  scale. 

The  asymmetry  produced  by  the  anisotropic  Knight  Shift  is 
clearly  seen. 
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With  the  line  being  narrowed  by  the  rotation,  it 
was  an  excellent  opportunity  to  accurately  measure  the  room- 
temperature  isotropic  Knight  Shift.  Previous  measurements 
have  been  complicated  by  the  anisotropic  Knight  Shift  (Borsa 
and  Barnes,  1964)  or  skin-depth  problems  (Jones  and  Williams, 
1964) .  The  applied  field,  corresponding  to  a  metallic  tin 
resonance  at  6  MHz,  was  determined  by  measuring  the  resonant 

frequency  of  an  undoped  saturated  aqueous  solution  of 

7  117  119 

LiCl.  Knowledge  of  the  ratios  of  the  Li  ,  Sn  and  Sn 

magnetogyric  ratios  lead  to  a  value  for  the  Sn  Knight  Shifts . 

117  119 

The  values  for  the  Sn  and  Sn  magnetic  moments  from  the 
literature  (Varian  Associates  NMR  Wall  Chart,  5th  ed. ,  1965)  gave 
values  for  the  Knight  Shifts  of  0.733%  and  0.747%  respec¬ 
tively.  This  difference  corresponds  to  a  hyperfine  structure 

-2 

anomaly  A  of  1.6  x  10  for  the  tin  nuclei,  in  contrast  to 

-4 

the  expected  value  of  the  order  of  10  (Eisenger  and 

117 

Jaccarino,  1958) .  Redetermination  of  the  ratio  of  the  Sn 
119 

and  Sn  magnetic  moments  removed  this  discrepancy.  The 
ratio,  determined  using  a  M/3  aqueous  solution  of  FeCl^ 
saturated  with  SnC^^  was  found  to  be 

119/y^  117  =  1.046535  +  0.000003. 

'  Sn  ' Sn  — 

Using  the  value 

YLi7/YSn119  =  1.0428555 

quoted  by  Borsa  and  Barnes  (1964) ,  the  isotropic  Knight  Shifts 


*  \ 
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117  119 

for  the  Sn  and  Sn  were  found  to  be  0.745  +  0.001% 
and  0.747  +  0.001%  respectively,  the  difference  being  attri¬ 
buted  to  experimental  error.  It  was  not  necessary  to  correct 
the  measured  resonance  frequencies  for  shifts  due  to  x'  since  the 
symmetry  of  the  line  indicated  a  negligible  contribution. 

To  determine  the  value  of  the  spin-spin  coupling 
constant  J  it  was  necessary  to  compare  the  calculated  spectra 
with  those  found  experimentally.  A  smooth  curve  is  drawn 
through  each  experimental  curve  and  is  fitted  against  the 
calculated  spectrum.  To  facilitate  the  drawing  of  this  curve, 
advantage  was  taken  of  the  symmetry  of  the  line  in  the  follow¬ 
ing  way.  Each  resonance  was  read  out  of  the  Fabri-Tek 
memory  twice.  One  graph  was  turned  relative  to  the  other, 
and  the  smooth  curve  drawn  through  the  superimposed  curves . 

This  has  the  effect  of  'reducing'  the  noise 

by  a  factor  of  /2.  The  experimental  curves  are  taken  in  the 
form  of  derivatives  and  for  this  reason  it  is  necessary  to 
calculate  the  theoretical  spectrum  in  the  same  form.  The 
computer  programme  for  calculating  the  spectrum  is  given  in 
appendix  III.  The  final  theoretical  curve  is  'assembled'  by 
weighting  the  spectra  from  the  various  configurations  of 
nuclei  according  to  table  I.  Corrections  were  made  to  these 
values  to  take  account  of  the  20%  of  the  nuclei  whose  spectra 
was  not  directly  calculated.  The  final  probabilities  are 
given  in  table  V.  Before  this  spectrum  is  compared  to  the 
experimental  curve,  it  is  broadened  to  allow  for  the  effect 
of  the  finite  field  modulation  and  rf  observing  power  (see 


Table  V 


Final  probabilities  of  the  various  nuclear  configurations 
used  to  calculate  the  theoretical  lineshape.  These  were 
used  instead  of  those  given  in  Table  I. 
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appendix  IV) .  The  spectrum  is  then  compared  to  the  experi¬ 
mental  curve  by  a  least  squares  fit. 

The  theoretical  curves  are  determined  by  the  para¬ 
meters  J^,  ,  J^,  B',  AW'  and  T^,  defined  in  the  first  part 

of  this  chapter.  It  is  also  necessary  to  allow  for  the 
differences  in  the  spin-lattice  relaxation  times  of  the  two 
isotopes.  This  may  be  appreciated  by  noting  that  the  T^ 
of  the  observed  nucleus  produces  the  lifetime-broadening 
effect,  whilst  the  T^  of  the  unlike  nuclei  produces  the 
collapse  of  the  satellite  lines  and  a  resultant  narrowing. 

It  is  difficult  to  determine  the  appropriate  values  of  the 
'broadening  factor'  B'  or  the  extra *f lip'  rate  AW'  and  these 
were  therefore  taken  as  unknowns.  It  was  found  that  taking 
reasonable  values  of  AW'  (-  0.05  W)  produced  little  change  in 
the  theoretical  lineshape.  This  result  could  be  anticipated 
by  noting  that  AW'  must  add  onto  the  effect  of  the  nuclear 
relaxation  processes.  This  parameter  was  therefore  set  as  zero. 
The  same  comment  cannot  be  made  for  the  factor  B'  which  has 

a  more  pronounced  effect  on  the  lineshape. 

117 

Of  the  two  isotopes,  Sn  should  lead  to  a  better 

value  of  the  spin-spin  coupling  constant.  This  is  because 

the  lifetime  broadening  is  smaller  (  (T^)Sn117  >  (T^)Sn119) 

117 

whilst  the  effect  of  the  J  coupling  is  larger.  (A  Sn 

nucleus  has  on  the  average  more  unlike  neighbours  than  a 
119 

Sn  nucleus  in  tin  metal.) 

In  the  calculation  of  the  theoretical  curves  it  is 
necessary  to  know  or  assume  the  ratios  between  the  first, 
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second  and  third  nearest  neighbour  couplings.  Whilst  it 
would  be  desirable  to  treat  all  these  values  as  unrelated 
parameters  to  be  determined  separately,  this  approach  leads 
to  computational  problems  and,  more  serious,  is  not  justified 
by  the  present  s ignal-to-noise  ratio.  Two  models  were  used. 

In  the  first  model,  it  is  assumed  that  only  the  second- 
nearest  neighbours  interact,  and  there  is  therefore  only  one 
parameter  to  be  determined  (J^,  J^,  B'  and  AW  are  taken 
as  zero) .  This  model  was  considered  in  view  of  the  suggestion 
made  by  Sharma  et  al.  (1969)  that  the  second  nearest  neigh¬ 
bour  interaction  is  dominant  on  account  of  the  shape  of  the 
tin  Fermi  surface.  This  model  gives  a  value  of  3  kHz  but 
the  fit  is  poor.  The  theoretical  curve  is  not  as  broad  as 
the  experimental  one,  peaking  too  early  and  fading  too 
quickly  in  the  wings.  In  the  second  model,  it  is  assumed 
that  the  couplings  have  the  form  k  F(x)  where  F(x)  is  the 
Ruderman-Kittel  relationship  (see  chapter  3)  and  k  is  a 
constant.  The  ratio  of  the  J  values  are  therefore 

J1  :  J2  :  J3  ::  1.00  :  0.440  :  -0.631, 
and  the  extra  broadening  factor  B'  has  the  form 

B'  =  BJ(2tt2J2)/W. 

The  undetermined  constant  BJ  is  approximately  0.02  when  one 
takes  the  typical  values  of  the  tin  parameters  given  in  the 
first  part  of  this  chapter.  Using  the  second  model  it  was 
found  that  the  first  nearest  neighbour  interaction  in 
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Sn  had  the  values  of  -2.21,  -2.04,  -1.88,  -1.76  and  -1.65 

(+  0.09)  kHz  when  BJ  was  0.00,  0.01,  0.02,  0.03  and  0.04 

respectively.  The  negative  sign  of  the  coupling  is  inherent 

119 

in  the  Ruderman-Kittel  model.  Values  from  the  Sn  '  resonance 
are  -2.40,  -2.09  and  -1.94  (+  0.09)  kHz  when  BJ  was  0.00, 

0.01  and  0.02.  As  can  be  seen  from  figure  (4.6)  the  discre¬ 
pancy  between  the  experimental  and  theoretical  resonances 
is  of  the  order  of  the  reading  error  in  determining  the  values 
of  the  smoothed  experimental  curves.  Figure  (4.7)  shows 
the  theoretical  curves  superimposed  on  an  experimental 
resonance . 

It  would  be  possible  to  calculate  the  theoretical 
curves  corresponding  to  J  values  obtained  from  models 
other  than  those  considered.  However,  consideration  of 
more  detailed  models  is  not  justified  by  the  present  signal- 
to-noise  ratio. 
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Figure  (4.6) 

A  graph  of  the  normalized  rms  discrepancy  between  experiment 

117 

and  theory  for  a  number  of  Sn  resonance  plots. 

The  parameter  BJ  was  taken  as  0.02.  The  normalization  is 
obtained  by  taking  the  peak-to-peak  height  of  the  experi¬ 
mental  resonances  equal  to  2.  The  dashed-line  gives  an 
indication  of  the  error  to  be  expected  if  the  smoothed 
experimental  curve  differed  everywhere  from  the  calculated 
spectra  by  an  amount  equal  to  a  reading  error  of  0.25  mm. 


:* 


. 


X 


■*’ 


NORMALIZED  R.M.S.  DISCREPANCY  ( arbitrary  units) 


J}  (kHz) 


83 


Figure  (4.7) 

117 

Experimental  Sn  resonance  obtained  by  making  1250  runs 
over  a  period  of  29  hrs .  The  sample  was  spun  at  5  kHz  and 
the  resonance  frequency  was  6  MHz.  The  frequency  scale  is 
549  Hz/cm.  The  smooth  curve  is  a  theoretical  fit  with 
=  1.75  kHz  and  BJ  =  0.02. 
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Chapter  5.  Discussion  and  Conclusion 

In  this  thesis  we  have  given  details  of  the  con¬ 
struction  of  a  high-speed  gas  turbine.  This  turbine  was 
capable  of  rotating  nmr  samples  at  the  'magic'  angle 
54°44»  within  a  Varian  crossed-coil  probe  at  speeds  in 
excess  of  5  kHz.  Rotation  at  the  magic  angle  at  speeds 
greater  than  the  equivalent  frequency  spread  of  the  reso¬ 
nance  linewidth  removes  the  effects  of  the  anisotropic 
Knight  Shift,  dipolar  interactions  and  the  pseudo-dipolar 
part  of  the  electron-coupled  nuclear-nuclear  interactions. 

Measurements  of  the  narrowed  resonances  gave  accurate  values 

117 

of  the  room- temperature  isotropic  Knight  Shifts  of  the  Sn 
119 

and  Sn  isotopes.  In  addition,  the  removal  of  the  aniso¬ 
tropic  interactions  enabled  the  pseudo-exchange  coupling  in 
tin  metal  to  be  measured.  The  spectrum  of  the  rotating 
sample  is  determined  by  a  'liquid-like'  Hamiltonian,  and 

enables  the  use  of  techniques  usually  more  appropriate  to 

117 

high-resolution  nmr  to  compute  the  lmeshape  of  the  Sn 
and  Sn'*"'^  resonances.  It  was  found  that  the  lineshape 
is  partly  determined  by  the  effect  of  the  finite  life¬ 
time  of  the  nuclear  spin-states  caused  by  the  spin- 
lattice  relaxation  time  in  tin  at  room  temperature.  It 
would  be  interesting  to  repeat  the  measurements  and  remove 
the  broadening  effect  of  the  short  T^  by  spinning  the  sample 
at  temperatures  near  liquid  nitrogen  or  lower.  This  could 
be  achieved  by  driving  the  turbine  using  cooled  helium  gas. 


~  \ 

j  • 


85 


However,  it  would  not  be  possible,  because  of  space  limita¬ 
tions  inside  the  Varian  probe,  to  convert  the  present  system 
to  such  an  operation.  Two  turbine  systems,  capable  of 
operating  at  low  temperatures,  have  been  described  in  the 
literature  by  Beams  (1971)  and  Potter  (1971).  Potter's  tur¬ 
bine  may  even  be  propelled  by  liquid  Helium  to  give  operation 
at  temperatures  near  1.2°K,  although  with  a  limited  speed 
range.  Both  these  turbines  differ  in  construction  from  the 
turbine  described  in  this  thesis.  The  propelling  and 
supporting  gas  streams  are  separated  from  each  other.  This 
offers  the  advantage  that  at  high  flow  rates  of  the  propelling 
gas  stream,  the  turbine  does  not  shift  its  position,  so  that 
the  flutes  and  gas  stream  remain  fixed  with  respect  to  each 
other.  This  provides  more  effecient  use  of  the  propelling 
gas  stream  and,  if  the  turbine  remains  balanced,  higher 
rotation  rates.  If  construction  of  low  temperature  turbines 
was  attempted  it  would  be  advisable  to  prepare  the  sample 
container  as  large  as  possible  to  obtain  the  maximum  signal- 
to-noise  ratio.  The  problem  of  poor  s ignal-to-noise  ratio 
was  one  that  placed  the  major  limitation  on  the  accurate 
determination  of  the  pseudo-exchange  coupling  in  our 
measurements. 

The  measurements  of  the  Knight  Shifts  made  when  the 

metallic  tin  resonance  was  narrowed  by  rotation  indicated, 

at  first,  that  there  was  a  significant  difference  in  the 

117  119 

values  for  the  Sn  and  Sn  isotopes.  However,  a  re- 

117  119 

determination  of  the  ratio  of  the  Sn  and  Sn  magneto- 
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gyric  ratios  removed  this  difference.  Our  measured  value 
of  the  ratio  was 

Y119/Y117  =  1.046535  +  0.000003, 

in  contrast  to  previous  values  of  1.04639  (Lindgren,  1965) 

and  1.04648  (Varian  NMR  Table,  fifth  edition,  1965).  The 

resonance  frequencies  of  the  rotated  tin  samples  were 

measured  relative  to  both  an  undoped  LiCl  solution  and  an  undoped 

SnCl2  solution.  Using  the  LiCl  reference  sample  and  the 

ratio 

V  7  /v  11p  =  1.0428555 

Li  (LiCl)  Sn  (SnCl2) 

quoted  by  Borsa  and  Barnes  (1964)  leads  to  values  of  the  room- 

117  119 

temperature  isotropic  Knight  Shifts  for  Sn  and  Sn  '  of 
0.745  +  0.001%  and  0.747  +  0.001%  respectively.  This  is  in 
good  agreement  with  the  room-temperature  value  of 

119 

0.750  +_  0.002%  obtained  by  Borsa  and  Barnes  for  the  Sn 

isotope.  However^ we  have  redetermined  the  ratio 

v  /v  q  using  undoped  aqueous  solutions  of  LiCl  and 
Li  Sn 

SnCl2  and  found  the  value 

v  7/v  119  =  1.04264  +  0.00003. 

Li '  Sn 

which  gives  the  Knight  Shifts  of  0.724+0.003%  and  0.726+0.003% 
which  are  in  closer  agreement  with  the  value  of  0.735+0.002% 
obtained  at  room  temperature  by  Sharma  et  al  (1969) . 

Dr.  Williams  (private  communication)  has  indicated  that  their 
measurements  were  made  with  direct  reference  to  a 
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paramagnetically-doped  aqueous  solution  of  SnCl2,  without 

correcting  for  the  effects  of  paramagnetic  susceptibility. 

7 

It  appears  from  our  measurements  that  the  ratio  of  Li  and 
119 

Sn  resonance  frequencies  taken  by  Borsa  and  Barnes  is 
incorrect  to  about  2  parts  in  10^.  Using  our  corrected 
ratio  the  Knight  Shift  measurements  of  Borsa  and  Barnes 
come  in  close  agreement  with  both  our  measurements  and  those 
of  Sharma  et  al. 

The  synthesis  of  the  theoretical  lineshape.,  correspond¬ 
ing  to  a  'liquid-like'  Hamiltonian^,  by  techniques  of  high 
resolution  nmr  spectroscopy  was  made  possible  by  the  small 
abundance  of  Sn^^  and  Sn^^.  The  analysis  showed  that  distant 
unlike  nuclei  had  a  significant  effect  on  the  resonance 
lineshape.  Taking  a  value  of  0.02,  a  reasonable  value  of 
the  parameter  BJ  describing  the  effect  of  distant  unlike 
neighbours,  we  have  found  experimentally  that  the  nearest- 
neighbour  pseudo-exchange  coupling  is  -1.89  ±  0.09  kHz. 

This  can  be  compared  with  the  values  of  2.5  kHz  (Karimov  and 
Schegolev,  1961)  and  2.0  +  0.5  kHz  (McLachlan,  1968)  obtained 
using  second  moment  techniques.  The  latter  values  are  rather 
high  due  possibly  to  the  neglect  of  the  far  neighbour  inter¬ 
actions.  The  broadening  effect  of  these  interactions  means 
that  a  smaller  value  of  pseudo-exchange  coupling  is  needed  to 
account  for  the  width  of  the  resonance  line.  Alloul  and 
Deltour  (1969)  obtained  a  value  of  -4.1  kHz  using  a  pulse 
method.  This  is  of  great  interest  as  it  should  in  principle 
lead  to  measurements  of  greater  accuracy  than  the  second  moment 
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method.  The  technique  in  the  pulse  method  is  to  examine  the 
spin-echo  envelopes  of  the  tin  nuclei  in  tin  alloys .  The 
presence  of  the  impurity  atoms  leads  to  a  spatial  distri¬ 
bution  of  tin  Knight  Shifts  related  to  the  oscillation  of 
the  charge  density  of  the  conduction  electrons  (Froidevaux 
and  Weger,  1964).  The  spatial  distribution  of  Knight  Shifts, 
provided  they  are  large  enough,  makes  neighbouring  nuclei 
dissimilar,  even  if  they  are  the  same  tin  isotope.  The 
coupling  between  unlike  nuclei  leads  to  a  modulation  of 

the  spin-echo  envelope.  This  was  clearly  shown  in  the  en- 

195 

velopes  obtained  from  Pt  '  in  platinum  alloys  (Froidevaux 
and  Weger,  1964).  In  contrast  to  platinum,  tin  has  two 
magnetically  active  isotopes ,  so  that  there  should  be  a 
small  modulation  in  the  spin-echo  envelope  of  pure  tin  metal. 
There  is  some  evidence  of  such  a  modulation  in  the  data  shown 
by  Alloul  and  Deltour  (1969)  though  they  do  not  draw  attention 
to  it.  The  modulation  amplitude  in  tin  increases  on  the 
addition  of  the  impurity  but,  even  so,  the  modulation  is 
not  as  dramatic  as  it  is  in  the  platinum  alloys.  This  is 
partly  due  to  the  difficulty  in  producing  stable  solid  solu¬ 
tions  with  sufficient  impurity  concentrations  to  provide  a 
large  distribution  of  Knight  Shifts.  In  platinum  alloys 
the  spin-echo  envelope  modulation  shows  fine  structure 
(Froidevaux  and  Weger,  1964)  but  in  the  tin  alloys  the 
modulation  dies  away  so  fast  that  only  the  fastest  modula¬ 
tion  period  will  be  seen.  Alloul  and  Deltour  have  attributed 
the  tin  modulation  to  nearest-neighbour  interactions. 
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Sharma  et  al.  (1969)  have  suggested  that  it  in  fact  arises 
from  a  dominant  second  nearest-neighbour  interaction.  They 
base  this  upon  the  presence  of  a  flat  region  on  the  tin  Fermi 
surface  in  the  [D01]  direction  which  could  lead  to  signifi¬ 
cant  deviations  from  the  Ruderman  -Kittel  interaction  (Roth, 
Zeiger  and  Kaplan,  1966) .  We  have  shown  in  this  thesis  that 
the  assumption  of  a  second  nearest-neighbour  interaction  far 
stronger  than  the  other  interactions  does  not  give  good 
agreement  with  the  data  obtained  with  the  resonances  from 
the  spinning  tin  samples.  It  is  therefore  necess¬ 
ary  to  offer  an  explanation  for  the  discrepancy  of  a  factor 
of  2  between  Alloul  and  Deltour  and  ourselves.  Froidevaux 
and  Weger  (1964)  have  shown  that  the  modulation  of  the  spin- 
echo  envelope  is  determined  by 

» 

E(2t)  =  £  A  cos  (2ttJt)  , 

r=0  r 

where  A  is  the  probability  of  having  r  nearest  neighbour 
magnetic  nuclei.  The  probability  of  having  two  magnetic 
nearest  neighbours  is  close  to  that  of  having  one  magnetic 
nearest  neighbour.  In  the  case  of  two  neighbours,  the  enve¬ 
lope  will  have  twice  the  modulation  period  than  where  there 
is  only  one  nearest  neighbour.  Because  of  the  rapid  decay 
of  the  modulation  and  the  poor  s ignal-to-noise  ratio,  it  is 
probable  that  only  the  faster  modulation  period  was  observed, 
giving  rise  to  a  period  determined  by  2J  rather  than  J. 

This  interpretation  of  Alloul  and  Deltour' s  data  would 
bring  their  value  of  J  into  closer  agreement  with  ours. 
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However,  it  must  be  pointed  out  that  the  analysis  of  Alloul 
and  Deltour' s  data  is  complicated  by  the  fact  that  in  an 
unspun  sample  it  is  necessary  to  consider  the  pseudo-dipolar 
and  dipolar  interactions.  These  have  the  effect  of  modify¬ 
ing  the  frequency  of  the  modulation.  It  is  also  interesting 
to  note  that  Froidevaux  and  Weger  (1964)  mention  that  their 
results  on  platinum  alloys  gave  results  which  were  a  factor 
of  two  greater  than  those  determined  by  steady-state  measure¬ 
ments.  A  second  possible  explanation  of  the  discrepancy  may 
be  offered  if  it  is  assumed  that  the  modulation  frequency  has 
been  correctly  assigned  by  Alloul  and  Deltour.  The  frequency 
is  determined  by  both  J^,  the  nearest  neighbour  scalar 
coupling  constant,  and  ,  the  nearest  neighbour  'dipolar 
and  pseudo-dipolar'  coupling  constant.  To  determine  these 
quantities  separately  it  is  necessary  to  have  an  additional 
relationship  between  them.  Alloul  and  Detour  calculated 
that  M the  theoretical  second  moment  of  the  resonance  line- 
shape  ,  is  given  by 

2  2 

M2  =  0 . 40B1  +  0.133^  . 

As  we  have  previously  mentioned,  the  experimental  determination 

of  the  second  moment  is  rather  difficult.  Alloul  and  Deltour 's 

2 

value  of  2.5±0.3  (kHz)  can  be  compared  with  the  value 

2 

1. 2±0.2  (kHz)  obtained  by  Karimov  and  Shchegolev  (1961). 

Using  the  modulation  period  determined  by  Alloul  and  Deltour 
and  our  value  of  =  -1.9  kHz  we  find  that  B^  = 


-1.8  kHz . 
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These  values  lead  to  a  second  moment  of  1.8  (kHz)  ,  which  is 
between  the  two  experimental  values.  Taking  this  value  of 
B^,  it  follows  that  the  pseudo-dipolar  coupling  between 
nearest  neighbours  is  -2.3  kHz.  It  must  be  pointed  out  that 
the  negative  sign  of  J-^  is  inherent  in  the  Ruderman-Kittel 
model  and  cannot  be  determined  from  the  calculated  spectrum. 

If  we  take  the  positive  sign  for  ,  then  this  leads  to  a 
pseudo-dipolar  coupling  constant  of  1.3  kHz. 

It  should  be  practical  to  repeat  Alloul  and  Deltour’s 
measurements  in  alloys  for  cases  where  the  contribution  from 
has  been  removed  by  macroscopic  rotation  or  by  a  multiple 
pulse  sequence.  This  would  lead  to  a  spin-echo  evelope 
determined  by  only  and  it  would  therefore  be  hoped  that 
a  more  accurate  value  of  the  pseudo-exchange  coupling  would 
be  obtained. 

In  conclusion,  we  suggest  that  it  would  prove  useful 
to  attempt  to  obtain  confirmation  of  the  J  values  obtained 
in  this  thesis  by  applying  the  second  moment  technique  to 
the  resonances  of  the  spinning  samples.  It  would,  however, 
be  necessary  to  correct  the  lineshape  for  the  broadening 
effect  of  the  T^  of  the  nuclei  being  observed.  It  is  not 
necessary  to  remove  the  effect  of  T^  on  the  neighbouring 
nuclei,  as  this  interaction  produces  relaxation  narrowing, 
which,  in  a  similar  fashion  to  the  motional  narrowing,  leaves 
the  second  moment  unchanged.  However,  the  narrowing  will  tend 
to  produce  a  lineshape  that  is  somewhat  Lorentzian,  with  the 
associated  difficulty  in  making  accurate  determination  of  the 


second  moment. 
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Appendix  I:  Calculation  of  the  AA'  part  of  the 

AA'X  spectrum. 

The  r^^'s  satisfy  the  simultaneous  equations  given 
in  table  IV.  The  solution  of  the  equations  is  not  diffi¬ 
cult  as  the  matrix  A  has  a  number  of  special  properties . 

It  is  symmetric,  and  both  the  imaginary  parts,  and  the  only 
quantity  which  varies  at  each  point  of  the  spectrum,  lie 
along  the  diagonal.  The  equations  for  r  may  be  written 

A  .  r  =  H  i  Z  P (VK) _1  (1) 

=  xo  o—  ' 

where  A  and_P_  are  given  in  table  IV.  The  nuclear  suscep¬ 
tibility  x(^)  has  ^e  form 

H  X(w)  =  r  .  P  yh  (2) 

A  U 

Let  r  =  C  +  i  D,  so  that 

H  x"(a))=“D*I)Yfr*  (3) 

xo  A  —  — 

Splitting  matrix  A  into  real  and  imaginary  parts  we  have 
£  =  F  +  i  G  (4) 

so  that  D  satisfies 

(£  +  Q  F~1  g)  D  =  P.  (5) 

It  is  necessary  to  solve  for  D  for  each  value  of  oj.  This 

is  easy  since  F  is  not  a  function  of  w  and  G  is  a 

diagonalmatrix  and  can  be  expressed  as 


. 


X 
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G  =  (k  -  w)  I 


(6) 


where  I  is  the  unit  matrix 


(7) 


and 


03  =  0)1  . 


For  computational  ease  it  is  better  to  keep  both  A  and  A' 
parts  of  the  spectrum  (which  gives  a  symmetric  solution) 
rather  than  just  the  A  line.  This  is  justified  as  it  will 
be  necessary  later  to  take  account  of  the  A'  resonance. 

To  avoid  including  the  A1  nuclei  twice,  the  relative  weight 
given  to  the  AA1  spectrum  in  determining  the  main  spectrum 
must  be  decreased  by  a  factor  of  two  over  the 
probability  for  finding  an  A  nucleus  at  the  centre  position 
given  in  table  I,  column  3. 


Experimentally,  it  is  that  is  of 


interest.  It  is  then  necessary  to  solve  for  9D/9a)  which 
satisfies  the  equation 


(F  +  G  F  1  g)D  =  E  1  §)  )  D  . 


(9) 


It  is  necessary,  first,  to  solve  equation  (5)  for  D  and  then 
substitute  these  values  into  equation  (9) . 


. 

w*  X 

- 
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•  *h 

Appendix  II:  Proof  that  the  antisymmetric  part  J  of  the 

indirect  nuclear  spin-spin  tensor  in  white 

tin  is  zero. 

To  have  a  spectrum  that  is  determined  by  the 
'liquid-like'  Hamiltonian,  it  is  necessary  that  the  anti- 
symmetric  part  JT  of  the  electron-coupled  nuclear  spin-spin 
interaction  tensor  be  zero.  Since  the  point  group  symmetry 
for  white  tin  is  4  2  m  (Pearson,  1967) ,  it  is  easy  to  show 
that  JT  must  be  zero. 

Consider  a  general  tensor  T  referred  to  the  usual 
crystallographic  axes  of  tin.  Applying  4  symmetry  requires 
that 


T  = 


T  T 

11  12 


-T  T 

±12  22 


0 


0 


0 


0 


33-J 


whilst  m  symmetry  further  requires  that 


T  = 


11 


0 


0 


0 


22 


0 


0 


0 


33 


This  implies  that  any  second-rank  tensor  describing  a 
physical  property  of  white  tin  must  be  symmetric  and  J 


f 


must  be  zero. 


* 


: 
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APPENDIX  III  :  Program  for  calculating  the  theoretical  lineshape. 


CURVE  {  1  ,  50  ) 

1  C  CALCULATION  OF  THE  LORENTZ-LIKE  LINE-SHAPE  UNDER  COEDITIOES  OF 

J  COUPLING , RELAXATION  HARROWING  AND  RELAXATION  BROADENING 
THE  PROGRAM  WILL  ACCOMADATE  BOTH  THE  CALCULATION  OF  THE  REAL  AND 
IMAGINARY  PARTS  OF  THE  NUCLEAR  SUSCEPTIBILITY 
THE  LINE SHAPES  ARE  GENERATED  IN  THE  DERIVATIVE  FORI 1 
ONLY  ONE  HALF  OF  THE  LINE-SHAPE  IS  GENERATED  AS  THE  RESONANCE 
IS  SYMMETRIC 


C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 
124 
25 
26 
27 
2  8 

2  9 

3  0 

b  i 

2 

3 

4 
i5 
6 
1 
8 
9 

It  0 
1 
2 
3 
•'■4 
|l  5 
6 
in 
8 
9 
0 

F  FILE 


C 
C 
C 

c 

904 


THE  PROGRAM  TAKES  INTO  ACCOUNT  THE  10  NEAREST  NEIGHBOUR 
LATTICE  SITES  IN  TIN  THESE  SITES  ARE  SPLIT  INTO  THREE  SHELLS 

OF  4,2  AND  4  SITES  RESPECTIVELY 

THE  PROGRAM  GENERATES  SPECTRA  TAKING  INTO  ACCOUNT  UP  TO  TWO 
OTHER  ACTIVE  NUCLEI  IN  THE  OUTER  SNELLS 

COMPLEX  CB 

DIMENSION  CB (1600), CK (1600) 

DIMENSION  G ( 1 9 )  , GG ( 1 S )  , X ( 1 0  C ) 

NUMBER  OF  THEORETICAL  CURVES  TO  BE  CALCULATED  IS  NRUN 
NUMBER  OF  POINTS  ON  THEORETICAL  CURVE  IS  HUM 

READ ( 5,904) NRUN , N UM 
FORMAT ( 212) 

PI- 3 . 141593 
//  UM1  -l'J  0M 
N  UM  2  =  N  UM*  1  9 
DO  950  NT-1 ,11  RUN 


C 
C 
C 

90  2 

C 
C 
C 
C 
C 

C 

c 
c 


VALUE  OF  FIRS 


T  NEAREST  NEIGHBOUR  J -COUPLING  IS  ASSUMED  TO  BE  RJ  HZ 


REA. £  (  5 ,9  02  )RJ 
FORMATi 7 FI  0.4) 

SJI-  2 *P1*RJ 

IF  CALCULATION  IS  TO  BE  PERFORMED  ON  THE  BASIS  OF  SECOND  NEAREST 
NEIGHBOUR  INTERACTIONS  ONLY  SET  N 2=2  FOR  IDENTIFICATION 
OTHERWISE  SET  N 2=0 

READ(  5 ,904)7/2 

INTERVAL  BETWEEN  THEORETICAL  POINTS  IS  SET  AS  PP  HZ 


1 

C 

c 

c 


REALMS  ,  902  )PP 
PQ-2 *PI *PP 
DO  1  1=1, NUN 
X( I) =( I-1)*PQ 

SPIN-LATTICE  RELAXATION  TIME  IS  T 1  SECS 
READ (5 ,9  02  )2'1 


. 


- 
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CURVE (51,92) 


51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 
87 

8 

89 

90 

91 

92 


J/=0. 5/2*1 


£7 

£7 

£7 

£7 


J/J  AMY  PP  USED  ON  TWO  COUNTS 

THE  FIRST  IS  TO  ACCOUNT  FOR  TEE  DIFFERENCES  IN  THE  SPIN-LATTICE 
RELAXATION  TIRES  BETWEEN  THE  TWO  ISOTOPES 
£7  AND  SECOND  TO  TAKE  ACCOUNT  OF  THE  EFFECT  OF  LIKE  NUCLEI  IN  THE 
£7  CON  TIN  I  UR 
C 

C  EXTRA  BROADENING  FACTOR  B-BJ*2*PI*PI*J1*J1/W  TO  TAKE  ACCOUNT  OF  THE 
£7  UNLIKE  NEIGHBOURS  IN  THE  CON  TIN  I  UR 

C 

READ(b  ,902) WJ , B J 
W 1  =  W* ( 1  +  WJ ) 

WW=W1*W1 

C 

C  READ  IN  ASSUMED  RATIOS  OF  THE  RUDE  Eli  AN  -KITTEL  J  VALUES 
C  THE  RATIOS  FOR  INTERACTIONS  BETWEEN  UNLIKE  NUCLEI  ARE  FED  IN  VIA  C 
£7  RATIOS  FOR  UNLIKE  NUCLEI  VIA  GG  GG(1)  TO  £7£7(4)  ARE  ZERO 

C 


READ(  5 ,902)(£7(J)  ,7  =  1  ,19) 

READ(  5  ,9C2)(£7£?(J)  ,  T=  1 ,19) 

CALCULATION  OF  AX  SPECTRUM  AND  A  SPECTRUM 

CALL  AX(W,W1  ,  BJ  ,  SJ I  ,CB  ,1!  UR ,  N  UM1  ,  WW  ,G,X) 

CALCULATION  OF  THE  X  PART  OF  THE  AA'X  SPECTRUM 

CALL  XA  A  X  (  W ,  W 1  ,  .5  J  ,  WJ  ,SJI  ,CB  ,  N  UR ,  N  £//'!,  WW ,  £7,£-  £7 , 1 ) 

CALCULATION  OF  THE  A  A  '  P/ii?Y  OF  THE  AA'X  SPECTRUM 

CALL  AAAX(  W  ,  W 1  ,  BJ  ,  S«7T  ,  CB  J1UR  ,  P£/A:T1  ,  WW ,  £7  ,  £7£7  ,Z  ,  f/V ) 
Z7(7  9  64  JJ-1,NUM2 
CK(JJ ) -AIRAGi CB(JJ) ) 


£7 
£7 
(7 

£7 
C 

/or 

O 

£7 
£7 
C 

962 
964 
£7 

C  THIS  IS  TO  READ  OUT  THE  CALCULATED  SPECTRA  WITH  IDENTIFING  PARAMETERS 
£7 

CALL  D URP ( i? J  ,  A7 2  , PP  ,  T 1  ,  !7J  ,  PP  ,  CK  , P PP ) 

950  CONTINUE 
STOP 
END 


JF  FILE 


. 


■ 


\ 
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T  CURVE 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 


93,111) 

SUBROUTINE  AX  (  W ,  W 1  ,  BJ  ,  SJI ,  CB  ,  A'  iW 

COMPLEX  CB  ,  [/I  2  ,  i/3  4 

DIMENSION  CB  (1600)  ,(7(19)  ,7(100) 

$  =  2*i7+07l+PO  *  00 I  *  SJ I  *0.5/  jy 
2)0  2  0  0  —  1  , 4 
GJ  =  0(0)  *SJ  I 

IFiJ.EQ.l ) G 0  TO  25 
10(0(0 )  .270. 0.  )00  TO  22 
25  00  21  JJ=1,NUM1 


N  UM 1 


WI/,G,X) 


l/12  =  -Q-Z(JJ)*(  0.  ,1  .  ) 

034= (01 2+00* (  0  .  ,  1  .  )/2)*(012-  GJ  *  (  0  .  ,1  .  )/2)-WW 

CB  (  00  +  (  0  -  1  )  *0  0//)  =  4  *  (  1  -  2  *  (  i/1  2  -  W )  *  i/1  2  /  i/3  4  )  *  1  0  0  0  0  /  i/3  4 

21  CONTINUE 
GO  TO  20 

22  00  23  00=1,0 OAf  1 

2  3  00(00+(0-1  )  *N  U'M  )  =  (  0  .  ,  0  .  ) 

20  CONTINUE 


RETURN 

END 


OF  FILE 


T  CURVE  {  112 

,139) 

112 

S UB ROUTINE  XAA X (  W ,  W 1  , BJ  ,  00  , SJI  ,CB  ,N UM , 

113 

COMPLEX  BX  ,  OS  ,  i/4 5  ,  i/3  6  ,  i/i  2  ,  i/7  8  , 000 , 0  ,  TOP 

114 

DIMENSION  CB  (1600)  ,  6:  (  1  9  )  ,  GG  (19)  ,  J  (  1  0  0  ) 

115 

00  40  J- 5,10 

116 

00=0(0) 

117 

000=00(0 ) 

116 

IHGJ.EQ.O.  .AND.CGJ  .EQ.Q.)GC  TO  21 

119 

00  41  00=1,0001 

120 

Zl=6'OX*  (OO  +  OOO  )/2 

121 

X2-SJI* (CJ-CGJ) / 2 

122 

Q  =  2*  17+  2  *  01  +00  *SJI*SJI*  0.5/  W 

123 

BX=-Q-( 0. , 1 . )*X(00) 

124 

U 4  5  =  ( 0 .  ,1.  ) *0 1 +BB 

125 

03  6= (  0  .  ,  -1 .  )*X1+BK 

126 

U  12  =  ( 0  .  ,1 . )*X2+BX 

1  27 

U 7  8  =  ( 0  .  ,  -1  .  )*a2+BB 

L  2  8 

000  =  01 2  *  03  6  *  04  5  *  07  8  - 4  *  NW*BK*BX 

129 

TOP  =  4 *  i/1  *00-  04  5  *  03  6  -  01  2  *  07  8 

L  3  0 

0  =  -  (  2 *TCP+  8  *00*  ( f/1  -00 )  -  2*00* TOP*  (  2*00*  ( 

.31 

C /DEN) /DEN 

.32 

00(00+(0-1 )  *0000  =0*10000 

.33 

41 

CONTINUE 

.34 

GO  TO  40 

.35 

21 

00  23  00=1,0001 

.36 

23 

00 ( 00+  ( 0-  1 )  *00/./)  =  (  0  .  ,  0  .  ) 

.37 

40 

CONTINUE 

38 

RETURN 

39 

END 

OF 

FILE 

*Wi,*BK) 


' 

.  . 


. 


. 


. 
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T  CURVE  { 140, 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150  912 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 
17  7 
17  8 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

1  89  8  2 

OF  FILE 


69 

70 

71 
83 


189  ) 

SUBROUTINE  AAAX(  7/ ,  i/1  ,BJtSJItCBt  N  UK ,  N  UM 1  ,  7/7/,  G  ,GG  ,  X  ,  7  <7 ) 

COMPLEX  CB 

DIMENSION  CB{  16  00)  ,<?(  19)  ,GG(  19)  8 ,8  )  tAI(  8 ,8  )  tTA(  8 ,8  )  ,TC(  8  ,8  )  , 

CTB(  8 , 8  )  ,IZ(  8 ,8  )  ,TZZ(  8 ,8  )  ,SX(  100)  ,I(  100  ) 

DIMENSION  FF ( 8 ) %EE ( 8 ) , DD ( 8 ) , DDL ( 8 ) 

£C>  9  0  <7  =  11 ,19 
GGJ=GG(J) *SJI 
GJ=SQRT(GU)*GU)*SJI*SJI ) 

IF (  CJ  .  ££  .  0  .  .717/1 . 6'£<7  .  .  0  .  )  GO  TO  21 

1=0 . 5*(oY/FI(  0 . 2 5 *GGJ *GCJ +  £</ *GJ  )  ) 

SPP=0. 5*GJ*W1/D 

SMM-SPP 

CPM-SPP 

CPP-0 . 2b*GGJ/D 
C MM- -CPP 
SPM- -CPP*W1 
CM-SQRT  (  0  .  5+1777//  2  ) 

SM-SQRT ( 0 . 5 -CRM / 2 ) 

CP^SQRTi  0. 5+FPP/2) 

SP=SQRT( 0 . 5 -CPP l 2 ) 

7y$P  =  2  *  7/  +  3  *  7/1  /  2  +  7/1  *  CPP*  CPF  /2+BJ  *SJI*SJ I*  0.5/  7/ 

CPP  =  CPP  *7/1 
C MM- -CPP 
DO  71  1=1,8 
DO  10  M- 1,8 
AI (  L  ,  /•' )  =  0  . 

>1  (  L  ,  77 )  =  0  . 

IF(L.EQ.M)  GO  TO  69 
<7(9  TO  7  0 
11(1  ,77)  =1  . 

/I  (1 ,77)  =  -WQP 
CONTINUE 
CONTINUE 
A( 1,2) -CPU 
1(1,3)  =0777/ 

>1(1,6)  -SPM 

A ( 1 , 7 ) = -SUM  / 

/U  2 ,4) =CPP 
/l  (  2  ,  5  )  =  -  CP// 

>1(2,8)  --SPP 
A ( 3 , 4 ) -CP M 
A  (  3,5) -SMM 
A(  3 ,8  )  =-5P/-7 
1(4,6) -3PP 
>1(4,7)  -SPM 
>1(5,6)  -CPM 
>1(5,7)  =  C  MM 
>1(6,8)  =  CPP 
>i(7,8)=CP// 

CONTINUE 


' 
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T  CURVE ( 190 

,  239  ) 

190 

DO  74  L-l , 8 

191 

DO  73  M-lyL 

192 

TA(L  ,  11 )  -A(L  , M ) 

193 

TA  (M ,L)  =A(L  ,27) 

194 

73 

A(M  ,L) =A(L  ,  M ) 

195 

74 

COLT I RUE 

196 

FF ( 1 ) =  -0 . 25 *( 2*GJ+GGJ) -D 

197 

FF( 2 ) =FF( 1 )+0. 5 *GGJ 

198 

FF( 3 ) =FF( 1)+GJ 

199 

FF( 4)=PP( 2)+GJ 

200 

FF( 5 ) =PP( 1 ) +2 *D 

201 

FF( 6 ) =FF( 2 ) +2*D 

202 

FF(  7  )  -FF(  3  )  +2*12 

203 

PF(  8  )  -FF (  4) +2*27 

204 

EE  ( 1 )  -Sll-CK 

205 

EE ( 2 ) -SP-CP 

206 

EE  (  3  )  =  -  (027+527) 

207 

PP(4)=-(0P+5P) 

208 

PP(5)=PP(3) 

209 

PP(  6  )=PP(  4) 

210 

PP ( 7 ) =PP ( 1 ) 

211 

PP(8)=PP(2) 

212 

C 

213 

c 

CALCULATION  OF  THE  NORMAL  RESONANCE  CURVE 

214 

n 

O 

215 

c 

216 

c 

GELG  IS  A  SUBROUTINE  FROL!  TEE  FORTRAN  SCIENTIFIC 

217 

c 

218 

CALL  GELGiAI ,TA ,8 ,8 ,1 .5-6,152?) 

219 

00  7  6  5  =  1,8 

220 

DO  75  27  =1,8 

221 

TB(L  ,27)  =i4  (L  ,Af)  +PP(  L  )  *FF(M)  *AI(L  ,27) 

222 

75 

TC(  L  , 27 )  =/i  J(  L  , 2/ )*(FF(L)  +FF ( 27 )  ) 

223 

76 

CONTINUE 

2  24 

t 

DO  80  JJ-ltNUMl 

225 

XXJJ=X(JJ) 

226 

DO  7  9  5  =  1 ,8 

227 

00(L) =PP(L) *10000 

00 

CM 

DO  78  27= 1,8 

CD 

CM 

TZ(5 , AO  -TB(L ,77) - XXJJ*TC(L9M ) +XXJJ*XXJJ*AI(L ,77 

o 

CO 

78 

TZZ(L ,  77  ) =TZ(L,M) 

231 

7  9 

CONTINUE 

’32 

CALL  GELGiDD  9TZ  ,8,1,1  .27-7  ,IF7?) 

33 

270  196  5  =  1  ,  8 

3  4 

DDD(L) =0 . 

35 

5(9  1  95  2:  =  1  ,8 

36 

195 

5275(5)  =555(5)  -  (  2  *XXJ  J*AI(L  ,  N )  -  TC{L  ,27)  )*55(27) 

37 

196 

CONTINUE 

38 

C 

39 

C 

CALCULATION  OF  THE  DERIVATIVE  RESONANCE  CURVE 

OF 

FILE 

SUBROUTINE 


PA  CKA C-E 


’ 


'  X 


100 


CURVED  2  40 , 251  ) 


40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

OF  FILE 


C 


30 


8  0 

21 

23 

90 


CALL  GELCiDDD ,TZZ ,8 ,1 ,1 . E -  6  ,  IFF ) 

SX  (  JJ)=  1  *  ( Z2Z2Z7  ( 1 )  *EE  (  1 )  +DDD  (  2  )  *5/27  (  2  )  +DDD  (  3  )  *FF  (  3  )  +DDD  (  4  )  *EE  (  4 
C)+DDD( 5 ) *FF( 5 ) +DDLH  6) *FF( 6 ) +DDD(1 ) *FF( 7 ) +FFF( 0 ) *FF( 8 ) ) 

CF(  JV  +  (J-1)  *777/17)  =  (  0.  ,1.  )*SX(JJ) 

CONTINUE 
CO  TO  90 
10  23  11  =  1  ,772/771 
CB(JJ+(J-1)*NUM)=(0. ,0. ) 

CONTINUE 

RETURN 


CURVE ( 252,265) 


52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 
6  4 
65 

OF  FILE 


902 

904 

910 


Z  OFF  0  UTINE  D  UMP  ( F  </  ,  77  2  ,  FP ,  ,'Z7 1  9WJ  ,BJ  ,  CK  ,  77  2/77 ) 

DIMENSION  CK (1600) 

l/RITE  (  6 , 9  0  2)  RJ 

WRITE  (  6 ,904)77  2 

WRITE (6,902 ) PP 

WRITE ( 6 , 9 1 0 ) FI 

WRITE (  6,902)  IV  , BJ 

N UM2 - N UM* 1 9 

WRITER  6 , 910  )  (  CK(JJ)  tJJ- 1  ,777/772  ) 

FORMA F  (  7 FI  0 . 4 ) 

FOF/F4F(  213) 

FORMA  T  (  7  F 1  0 . 4  ) 

RETURN 

END 


_ 
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Appendix  IV  :  Correction  for  saturation  and  modulation  effects. 


Modulation  and  saturational  broadening  differ  in 
their  effect  on  the  resonance  line-shape.  The  saturation 
broadening,  which  is  produced  by  the  action  of  the  rf  field 
on  the  population  distribution  of  the  nuclear  spin-states, 
directly  affects  the  sample  whilst  the  modulation  broadening 
is  totally  instrumental.  We  therefore  correct  the  theoretical 
resonances  to  allow  for  the  saturation  effects  and  then  modu¬ 
lation  broaden  the  resulting  line-shape. 

The  nuclear  absorption  x"  >  in  the  presence  of  a 
linear  polarized  rf  field  2H^  and  frequency  v,  is  determined 
by  (Bloembergen ,  Purcell  and  Pound  (BPP) ,  1948) 


XVv.H.p  =  X"(V,0)[1  +  'l'2H^T1g  (v)/2] 


(1) 


where  x" (v,0)  is  the  unsaturated  susceptibility,  determined 
by  a  normalized  lineshape  function  g(v).  The  BPP  theory  is 
known  to  be  incorrect  for  solids  where  the  spin-lattice  inter¬ 
actions  are  much  weaker  than  the  spin-spin  interactions 
(Redfield,  1955).  This  is  not  the  case  in  metallic  tin, 

and,  in  this  case, the  saturated  derivative  lineshape  f  .  (v), 

3X" (v,H1)  Sat 

proportional  to  - ^ - ,  is  given  by 


fsat(v)  =  Ag’  (v)  [1  +  pHl  g(v)//g(Vmax)]  2 

2  2 

where  A  is  the  gain  of  the  system,  p  is  y  Tg9(vmax)/2  anc^ 

v  is  the  frequency  of  one  of  the  extrema  of  g' (v)  provided 
max 

w  ,  the  frequency  of  the  magnetic  field  modulation  satisfies 


'•1* 
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a)  T,  <<  1.  We  may  therefore  use  the  calculated  derivative 
ml 

lineshape  f(v)  to  determine  the  theoretical  'saturated' 
lineshape 

o  v  Vmax  , 

f  sat  (V>  =  f(v)tl  +  PH^  /f(v')dv'{/  f(V)dv'}  L]  (3) 

o  o 

The  proportionality  constant  p  must  be  determined  from 
experiment.  From  equation  (2)  we  have  that 


[f  .  (v  )  ]  =  [f  (V  )  ]  [1  +  pH?]  2 

sat  max  exp  max  exp  1 


(4) 


so  that  a  plot  of  the  peak-to-peak  intensity  of  the  derivative 

2 

line  as  a  function  of  H,  gives  p  .  Allowance  must  be  made 

1  ^  exp 

for  the  intrinsic  increase  in  the  signal  intensity  with  H^. 

2 

In  the  experiment  on  tin,.  pH^  was  0.045  ,  implying  that  the 
correction  to  the  lineshape  was  small. 

The  technique  used  to  calculate  the  modulation  broad¬ 
ening  is  that  given  by  Flynn' and  Seymour  (1960) .  The  modula¬ 
tion  peak-to-peak  amplitude  is  taken  as  2a  (in  frequency  units) , 
so  that 


V+a  2  2  1/2 

f  sat (v)  =  N/  (a  "  (v-V)  )  7  fsat(v')dV 
v-a 


(5) 


where  N  is  a  normalization  constant. 


— 
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